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The Symmetric Vibrations of Aircraft 


R. W. TRAILL-NASH, Ph.D. 


(Depariment of Supply and Development, Aeronautical Research 
Laboratories, Melbourne) 


(Summary: The eigenvalue equations for symmetric vibration of a complete 
aircraft are derived in a very general form. The “lumped mass” approximation to the 
continuous mass distribution is used and sub-matrices are associated with properties 
of relatively simple branches of the system. The final eigenvalue equations are 
expressed in terms of these sub-matrices, so that in a numerical application the 
physical system, as such, is considered only in relation to the properties of the simple 
branches. It is assumed initially that the aircraft wing and tail have flexural axes 
of the conventional type, but it is shown in the Appendix that under certain 
conditions the treatment can be generalised to cover swept and cranked wing aircraft.) 


Introduction 


The calculation of natural frequencies and modes of vibration of a complete 
aircraft is a complicated and laborious process. It is complicated because of the 
nature of the dynamical system, and laborious because any representation of the 
system must involve a relatively large number of degrees of freedom. The labour 
of computation required in the determination of frequencies and modes of a 
particular aircraft must be accepted as unavoidable since, whatever the method of 
solution, the amount of computation is considerable. The derivation of the vibration 
equations is a task which need be performed but once in general terms to be 
applicable to the whole field of conventional aircraft. It is the purpose of this paper 
to develop the eigenvalue equations for symmetric vibrations of the widely used 
“lumped mass” approximate representation of the continuous system, in a form 
which is convenient for direct substitution of numerical values of fundamental 
physical quantities of the system. 


The eigenvalue equations are stated in matrix notation") and the derivation 
of these equations is set out in matrix form. This leads to concise statement 
throughout and is an aid to clarity. Considerable use is made of sub-matrices which 
are related to specific parts or branches of the system, and the equations finally 
derived are defined in terms of these matrices. This implies that in an application 
to a specific aircraft reference is made to the physical system only in the description 
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of properties of relatively simple branches and that a purely formal process of 
matrix multiplication and addition then yields the required equations. Further, 
the matrices involved in general refer to quite arbitrary numbers of degrees of 
freedom in the various branches of the system so that it is the form only and not the 
size (in terms of degrees of freedom) of the system which is specified. Initially it 
is assumed that wings and tail have flexural axes, but it is shown in the Appendix 
that this is not essential to the treatment. 


The present treatment of the problem is in some ways similar to that proposed 
by J. Morris and G. S. Green’. The same concept of “branch system” is used 
but whereas Morris and Green determine the frequencies and modes of the branch 
system under prescribed conditions of constraint (e.g. wing, encastré at the root) 
as a preliminary step in establishing the frequency equation for the complete system, 
this treatment proceeds directly to a formulation of the eigenvalue equations for the 
complete system. The main point of variance in the two methods lies in the manner 
in which the equilibrium conditions are utilised. Morris and Green use the 
equilibrium conditions to establish the frequency equation when the vibration modes 
and frequencies of the branch systems are known, while here these conditions lead 
to the establishment of a generalised inertia matrix. It is not practicable to devote 
the time and labour to a numerical comparison of the two methods. It is felt, 
however, that solution by matrix iteration’ of the equation presented is probably 
more economical in labour than the method proposed by Morris and Green. This 
is especially so in the calculation of a small number (say three) of lower order modes 
and frequencies, only for which, as pointed out later, can calculation be accurate. 
If, on the other hand, it were practicable to calculate a large number of modes and 
frequencies the relative expediency of the two methods would be more doubtful. 


The “lumped mass” approximation to the mass distribution associated with 
a mass-less continuous elastic system, which is adopted here and by Morris and 
Green, is a representation which is commonly used for this type of problem, and 
one which yields satisfactory comparison with experimental results for modes of 
lower order when the physical properties of the system are determined with sufficient 
accuracy. It is well known that agreement between calculated and experimentally 
determined natural frequencies and associated modes becomes less satisfactory with 
increase in order of the mode. As pointed out by A. R. Collar’, one source of error 
lies in the common practice of neglecting the contribution of shear distortion of beam 
webs when flexibilities of the structure are calculated; but this source can be 
eliminated in the calculation of influence coefficients. Further, if influence 
coefficients are measured directly from the aircraft this question is automatically 
eliminated—and in fact this source of error is covered by the provision that “ the 
physical property of the system be determined with sufficient accuracy.” 


Another more serious source of error lies in the phenomenon of shear lag, 
which source: ©) it is virtually impossible to eliminate. Shear lag introduces changes 
in the stiffness of the system for each mode, the effective stiffness depending on the 
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THE SYMMETRIC VIBRATIONS OF AIRCRAFT 


shape of the mode in question; so that even if the stiffness is determined directly by 
static test the effective stiffness in higher order modes is not correctly represented. 
Although calculated higher order natural frequencies and associated modes have 
little relation to those of the actual system, this is not to say that valuable data 
cannot be obtained by calculation and that recourse must be made to some form of 
experiment in all cases. For dynamical problems in which natural modes of 
vibration with associated generalised co-ordinates may be used to advantage, the 
displacement of the system can often be expressed adequately in terms of a small 
number of lower order modes. 


2. Representation of Aircraft 


It is assumed that the continuous mass distribution of the aircraft is replaced 
by a finite number of discrete masses at selected points on the elastic structure. The 
magnitudes and positions of these masses are selected to give as close an approxi- 
mation as possible to the continuous distribution. A knowledge of the types of 
displacement which will occur in the various parts of the system during vibration 
indicates the types of inertia forces which are of importance, and consequently to 
some extent prescribes the distribution. In symmetric vibration the fuselage 
displacements are translational with a small component of rotation about axes 
parallel to the wing flexural axis. The inertia effect in rotation is known to be small, 
so that the fuselage mass distribution may be represented by point masses on the 
centre-line. On the wings and tail the effect of rotary inertia is similarly small so 
that it is satisfactory to represent the spanwise distribution by masses with zero 
moment of inertia about chordwise axes. The wings and tail, however, do rotate 
in the torsional sense, and since the mass is distributed over a finite chord, inertia 
forces in the form of torques result from acceleration in rotation. It is necessary 
therefore to adjust the chordwise distribution at spanwise mass stations to take 
account of the torsional inertia forces, and this can be done in either of two ways. 
At each spanwise station two point masses may be located on the chord so that their 
centroid lies at the centroid of the chordwise distribution and their moment of inertia 
corresponds to that of the chordwise distribution. Alternatively, an appropriate 
moment of inertia may be assigned to one mass situated at each chordwise centroid. 
The latter alternative is adopted here as this not only effects an algebraic 
simplification in the case where a straight flexural axis is assumed to exist, but also 
provides a convenient set of co-ordinates for a consideration of aerodynamic forces 
when natural modes of vibration are used in aeroelastic applications. 


The elastic properties of the continuous system are retained unaltered, but since 
the inertia properties are reduced to a finite number of discrete values it is no longer 
necessary to express the former in continuous form. Influence coefficients, defining 
the displacements in various mass co-ordinates due to forces appropriate to a 
particular co-ordinate under prescribed conditions of constraint of the system, then 
serve to specify the elastic properties. 
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Figure 1 depicts a typical idealisation of a four-engined aircraft for vibration 
analysis. It is assumed initially that the wings and tail have straight flexural axes. 
The fuselage is replaced by a straight mass-less beam with point masses attached, 
while the wings and tail are represented by straight mass-less beams of known 
flexural and torsional properties. The wing and tail masses are offset from the 
elastic beams (which are located at the flexural axes) and with each mass is 
associated a polar moment of inertia. A further polar moment of inertia is 
associated with the mass at the tail centre-section to take account of the inertia 
moment arising from the tail fin in single fin aircraft. 


3. Notation and Co-ordinates 


Since the system is free-free, i.e. is reactionless and has no fixity in space, it is 
necessary to refer elastic displacements to some arbitrary datum within the system. 
It is convenient, as will become evident, to select this elastic datum at the wing- 
fuselage junction station and to measure elastic displacements relative to the tangent 
plane at this point. Elastically the system is then divided into three; and the 
fuselage forward of the reference station, the wing, and the part of the system aft 
of the reference station are to some extent treated separately. This requires 
differentiation in notation between fuselage properties fore and aft of this point. 
The subscript “f” refers to quantities relating to the fuselage and the qualifying 
subscripts “w” and “t” indicate the location relative to the reference station; for 
example, the subscript “ft” refers to fuselage properties on the tail side of the 
reference station, while “fw” refers to those forward of the reference station. 
Properties of the wing and tail are associated with the unqualified subscripts “w” 
and “t” respectively. 

The eigenvalue equations are derived for symmetric vibrations only so that it 
is not necessary to consider the complete system but merely one half which lies to 
one side of the plane of symmetry. The plane of symmetry is a vertical plane on the 
fuselage centre-line, and this implies that fuselage properties be halved throughout. 
It will be understood therefore that inertia and stiffness properties associated with 
the fuselage refer to half the corresponding values for the complete aircraft. 


Section 3.1. 
As indicated by the choice of subscripts, for the definition of properties the 
system can conveniently be considered in four branches. These are: 
(i) the wing, 
(ii) the fuselage forward of the reference station, 
(iii) the fuselage aft of the reference station, 
and (iv) the tail, 
and the following definitions are in accordance with this branch notion: 
(m,). | mass at i wing station 
(pi) moment of inertia at i wing station 
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Idealised representation of aircraft. 


(mje mass at i™ station of fuselage forward of reference station 
m, mass at reference station 

(mj: mass at i station of fuselage aft of reference station 
Px moment of inertia at fuselage-tail junction 
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(m,) mass at i™ tail station 
(p): moment of inertia at i™ tail station. 


Section 3.2. 


The system is referred to Cartesian co-ordinates, the x co-ordinate measured 
parallel to the fuselage centre-line and positive in the tail-to-nose sense, the z-axis 
downwards, and the y-axis completes a right-hand system. In the undistorted 
condition the aircraft representation is assumed coplanar and is positioned to lie in 
the xy-plane with reference station at the origin. 


Further constants relating to position in the x direction are then defined as 
follows : 

(rw distance of mass at i wing station from wing flexural axis 

(r;, | distance of mass at i™ tail station from tail flexural axis 
(in both cases measured positive in the positive x direction) 

(S;)iw | xX CO-ordinate of mass at i“ station of fuselage forward of reference 
station 

(si): | x co-ordinate of mass at i™ station of fuselage aft of reference station 

x co-ordinate of mass at i“ wing station 

(s;; | x co-ordinate of mass at i™ tail station. 


Section 3.3. 


Displacements of the system initially appear in the equations of motion in two 
forms—one as absolute displacements which are associated with inertia forces, 
and the other as relative displacements within the system which are associated with 
elastic forces. The absolute displacements are measured relative to the fixed 
co-ordinate axes and the notation is as follows* (see Fig. 2): 


(hi)~ absolute linear displacement of i“ wing mass 
(¢) absolute angular displacement of i” wing mass 
(hi absolute linear displacement of i mass of fuselage forward of 
reference station 
(h;) | absolute linear displacement of i“ mass of fuselage aft of reference 
station 
gr absolute angular displacement at fuselage-tail junction 
(h;) absolute linear displacement of i™ tail mass 
(v; | absolute angular displacement of i“ tail mass 
h, absolute linear displacement at reference station 
¢, absolute angular displacement at reference station. 


*In accordance with the right-hand system of axes, linear displacements are measured positive 
downwards and angular displacements positive in the nose-up sense. 
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WING FUSELAGE 


JUNCTION STATION 0 
x 
(s) 
Po Si 
(Sides G)= 9, tH) 
z 
Fig. 2. 
Displacement co-ordinates. 
Section 3.4. 


The displacements relevant to elastic forces are measured relative to the tangent 
plane at the reference station as shown in Fig. 2. These relative displacements are 
termed “elastic displacements ” and in the same order as the absolute displacements, 
the co-ordinates are defined as 

(u;)w elastic linear displacement of i“ wing mass 
(vi elastic angular displacement of i“ wing mass 


(uw elastic linear displacement of i mass on fuselage forward of 
reference station 


(ui) elastic linear displacement of i“ mass on fuselage aft of 
reference station 


vs elastic angular displacement at fuselage-tail junction 
(u; elastic linear displacement of i“ tail mass 


(v4. elastic angular displacement of i” tail mass. 


Section 3.5. 


The quantities so far defined refer to individual masses of the system and 
these are now grouped in matrix form to give a more concise notation. Where more 
than one matrix is associated with a particular branch of the system it is understood 
that the elements referring to specific stations occur in the same sequence in each 
matrix. The relevant matrices are defined in Table I. 
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The orders of these matrices depend on the representation of the system. In 
addition to those of the table the single term matrices pr, ¢, and ¥. are used and 
these are defined in Sections 3.1, 3.3 and 3.4 respectively. Further, certain 
“ operational” matrices are required and these are defined as 


I _ the unit matrix, 
La ccolumn matrix, each element of which is unity, 
H a square matrix, each element of which is unity. 


The orders of these matrices are nowhere specifically defined, and each may appear 
as matrices of different order in a single expression. This does not lead to confusion 
since they occur in association with matrices of known (although here, arbitrary) 
order, and hence the appropriate order of a particular operational matrix is implied 
by conditions of compatibility in the operations concerned. 


Section 3.6. 


It remains to define the elastic properties of the system. In the definitions 
which follow, square flexibility matrices are matrices of which the element in the i 
row and j” column represents the displacement in the i“ co-ordinate produced by 
a unit force in the j* co-ordinate under specified constraint conditions of the system. 
In this regard “displacement” and “force” may apply either to linear or angular 


TABLE I 
Elements | Elements 
Aw (maw hy (hiw 
P. (Diw | 
(mew | lew (hides 
Arn (mre hr (Ayr. 
A (m)) | h, (hi) 
Ry (re | Uw (uw 
R. (ri): Ve Ww 
Stw (Siew Utw (Us)tw 
Si: Ur (ui 
Sw (Si)w Ur (ui) 
Ve (Vi 


displacement and either linear force or moment respectively. 
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The conditions of 
| constraint are: for the wing and tail, encastré conditions at the root; for the fuselage 
fore and aft of the reference station, encastré conditions at the reference station. 


The matrices are defined as follows: 


flexibility matrix of fuselage forward of reference station 
flexibility matrix of fuselage aft of reference station 
linear dispiacement flexibility matrix on wing flexural axis 
linear displacement flexibility matrix on tail flexural axis 
angular displacement flexibility matrix of wing 

angular displacement flexibility matrix of tail 


column of linear displacements on fuselage aft of reference station 
due to unit moment at fuselage-tail junction 


(transposed of d) row of angular displacements at fuselage-tail 
junction due to unit linear forces at fuselage masses 

column of linear displacements on fuselage aft of reference station 
due to unit linear force at fuselage-tail junction (a column of Ey,~*) 


(transposed of e) row of linear displacements at fuselage-tail 
junction due to unit linear forces at fuselage masses (a row of E;,~*) 


linear displacement at fuselage-tail junction due to a unit linear 
force at this station (an element of E;,~*) 

linear displacement at fuselage-tail junction due to a unit moment 
at this station=angular displacement at fuselage-tail junction due 
to a unit linear force at this station 


angular displacement at fuselage-tail junction due to a unit moment 
at this station. 


4. Inertia Forces 


In the free-free vibration of a conservative system the only forces present are 
elastic and inertia forces. Elastic forces depend solely on the deformation of the 
system, while inertia forces are proportional to absolute accelerations. D’Alembert’s 
principle makes it possible to write down immediately the inertia force acting at 
any mass as 


force = — (mass) x (acceleration) for linear motion, 


moment= —(moment of inertia) x (angular acceleration) 


for angular motion. 
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(m Dee (p.), 


(Kis 
FLEXURAL 


(m,), (h ie 


2 
Fig. 3. 
Inertia forces at the i” wing station. 


In the case of simple harmonic motion at a frequency of w radians per second 


these become respectively : 
force = «? x (mass) x (displacement) 


moment = w? x (moment of inertia) x (angular displacement). 


and 
The elastic displacements in the various branches of the system resulting from 
inertia forces in vibration are now deduced in terms of these inertia forces. 
As a preliminary to the deduction of elastic displacements it is of advantage 
Fig. 3 depicts the 


to consider in detail the inertia forces on the wing and tail. 
conditions at a spanwise station of the wing. The linear inertia force at the i™ 


station is: 
(Fiw =? (miw (1) 
(2) 


and the moment about the flexural axis is 
=? { — (ri)w (mw (hiw } 


where F is measured positive in the positive z direction (downwards) and M is 


measured positive anti-clockwise. 
Similar expressions hold for forces at the i station of the tail. Further, at the 
tail root (fuselage-tail junction station) the linear force and moment from the tail 
10 
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Fig. 4. 
Inertia forces on fuselage aft of reference station. 


are derived by summation with respect to i. The force and moment from the tail 
at this point are 


F=w* (hi) (3) 

and M= { (pik Gi (dd 

or F=w? L’ A, h, ; (5) 


from Table I, where L’ is the transposed of L. 


5. Elastic Displacements 


The elastic displacements in various co-ordinates are as follows, where for 
convenience the fuselage-tail junction station is referred to as the point “C.” 


(i) Elastic linear displacements of fuselage aft of reference station 


Figure 4 illustrates the forces applied to this part of the fuselage. The elastic 
displacements consist of 

(a) displacements due to forces from fuselage masses, 

(b) displacements due to moment from moment of inertia at C, 
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(c) displacements due to linear force from tail at C, 
(d) displacements due to moment from tail at C. 


The displacement matrix with components in the order (a)—(d) becomes 


un = w? { E.,-* An hn t+dpu gn teL’ A, h,+dL’ (P, A; h.) } (7) 


(ii) Elastic angular displacement of fuselage-tail junction station 


The components of displacement are described as under (a)—(d) of (i). With 
components in this order the displacement becomes 


{ d’ An hy + 2D ornt+BL’ A. h.+aL’ (P, . (8) 


(iii) Elastic linear displacements of tail 


Displacements on the tail depend not only on displacements relative to the tail 
root but also on the linear and angular displacements at the root. The components 
of tail linear displacements are as follows 


(a) tail flexural axis displacements at spanwise mass stations relative to 
tail root, 


(b) displacements relative to root due to torsion about flexural axis, 


(c) displacement at C (this is included in (i) but requires to be specified 
separately. It consists of components as under (a)—(d) of (i)), 


(d) displacements due to angular displacement at C (the angular displacement 
at C is Vee of (ii) ). 


With components in the order (a)—(d) the displacement matrix then is 


u, =? [E.-* A. h.—R: (P, 
+L { An hi + on +8L’ A, h,+BL’ (P, A, h,) } 


(iv) Elastic angular displacements on tail 


The angular displacements on the tail similarly depend on the angular 
displacement at the root. The two components of displacement are 


(a) angular displacement relative to root, 


(b) angular displacement at C (v», of (ii)), 


12 


THE SYMMETRIC VIBRATIONS OF AIRCRAFT 


and with components in this order the displacements are 


[9,- vt —R, A, h.)+ 
+L {d An (10) 


(v) Elastic linear displacements of fuselage forward of reference station 


These displacements are due solely to linear inertia forces on the fuselage 
forward of the reference station. The elastic displacement matrix is then 


Ute =? Awe hen (11) 


(vi) Elastic linear displacements on wing 


These are of the same form as the displacements relative to the root on the tail. 
[he components are 


(a) wing flexural axis displacements 
(b) displacements due to torsion about flexural axis. 


The linear displacement matrix is 


(viii) Elastic angular displacements on wing 


Wing elastic angular displacements are measured relative to the root and are 
given by 


Ye= 07 : (13) 


The equations (7)—(13) define the elastic displacements of the system relative 
to the tangent plane at the reference station. More concisely they can be expressed 
in one matrix equation, and it is easily verified that they are equivalent to (14), in 
which broken lines indicate partitioning into sub-matrices. 


6. Equilibrium Conditions and Transformation of Co-ordinates 


The equation (14) contains two sets of displacement co-ordinates, the elastic 
displacement co-ordinates on the left hand side and the absolute displacement 
co-ordinates on the right hand side. These sets of co-ordinates are related through 
the equilibrium conditions of the system, i.e. the conditions of zero body displace- 
ments in translation and rotation; a consideration of these conditions yields the 
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transformation matrix for the two sets. This makes it possible to write (14) in 
terms of one set of co-ordinates to produce the required eigenvalue equations. 


The absolute displacements are expressible in terms of the elastic displacements 
and the displacements at the reference station. Thus, as shown in Fig. 2, for the 
i mass in a branch of the system, 


and «© «© «©  « (16) 


where the appropriate subscript is assigned to the terms in brackets. 


The complete set of relations (15) and (16) may be expressed as: 


h, 0 0 
ft Un St: L I L 
h, uy Si L I L 
wl °° Se L L 
or —] 0 | 


(17) 


In order to simplify the algebra certain matrices are now grouped together 
and re-defined. 


Let 


(18) 
Si: Ar Pr 
Stw Aw A Fe P 
Se Aw 


whence from (18), (17) becomes 


h, l l 
h 


hy Ure ca | 
hy Uv 
&& 
| L | L 
15 


R. W. TRAILL-NASH 


Section 6.1. 


For equilibrium the sum of all inertia forces in the z direction is zero and the 
moment about any point is zero. For the moment calculated about the reference 
station the equations of equilibrium become respectively 


[1 L’ m, h, 
Pile | 
and [Oo L’ 0 my h, 
S A h 0. (21) 


Substitution from (19) in (20) and (21) yields respectively 


L’ m, “| [0 [1 L’ 0} [m, 
A u A S Li+ 
P| | v P| | L 


+h, [1 L’ 0] 1] 
A I L \=0, 
a PIL 
(22) 
and [0 L’ Fim, “| fo 
S A 
= P| LY 
(0 L’ L] Pm, To 
A S Li+ 
+h, [0 L’ “| fim, k 
A I (23) 


Equations (22) and (23) are a set of two simultaneous equations in the variables 
h, and ¢, which are soluble in terms of the elastic displacements, i.e. the elements 
of u and vy. 


It is evident that the rows and columns corresponding to m, in the products 
involving u and ¥ can be omitted since they make no contribution to these terms. 
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The equations then respectively reduce to* 


(L’ ASL)<,—(m,+L’ AL)h,=(L’ 0) [ A 


and 


z 


(L’ SASL+L’ PL) =, —(L’ SAL)h, =(|L’ 
(25) 
The solution of these equations is 
LIT A u | A u ] 
and 
J PiLy @ 


where €=(L’ ASL)/y, 


o=(L’ SASL+L’ PL)/y, 


p=(m,+L’ AL)/y, 


y=(m,+L’ AL)(L’ SASL+ L’ PL)—(L’ ASLY , 


Section 6.2. 
Substitution for h, and <, from (26) and (27) in (19) then yields 


ap 


where the square matrix on the right hand side is the matrix of co-ordinate 
transformation. 


*The coefficients of the variables in (24) and (25) have the following physical interpretation : 
(m, + L’AL)= total mass of system, 
(L’ASL)= first moment of mass about reference station, 


and (L’SASL + L’PL)=second moment of mass about reference station. 
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It is now possible to express (14) in terms of the one set of co-ordinates { u,v } * 
and hence to obtain the required eigenvalue equations. In (14) the absolute 
displacement column { f,< } is pre-multiplied by the inertia matrix corresponding 
to'{h,¢}. On substitution of (29) this inertia matrix is modified by post- 


multiplication by the transformation matrix of (29), and the resulting product is the 


generalised inertia matrix corresponding to the co-ordinates { u,v }. 


If B denote this generalised inertia matrix, it is defined as 


[ B ] [ A ] < the transformation matrix 
L 


where use is made of the identities 
4], 


and H is defined in Section 3.5. 


The eigenvalue equations for symmetric vibration are then of the form: 


where D is the flexibility matrix in (14) and B is the generalised inertia matrix (30). 


7. Summary of Analysis 


The preceding analysis leads to a statement of the eigenvalue equations for 
symmetric vibrations in a concise form. For numerical work on a particular aircraft 
it is necessary merely to insert the appropriate values of the physical properties in the 


*Curly brackets thus { } denote a column matrix. 
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final equations and the problem is set up for solution without explicit consideration 
of dynamics. In this section the results of the analysis are summarised to provide 
ali the necessary information for the calculation of the eigenvalue equations. 


The equations are 


Vn 
uy 
where v=| 
Utw 
Yw 


and the terms with subscripts are defined in Table I and Section 3.4; D is the 
flexibility matrix (first square matrix on right hand side of (14)); B is the generalised 
inertia matrix (30), and in respect to (30), 


Ar Si: p 
ft 
A=| “4 , P=| PB 
Aww Stw P 
Aw S = 


w 


and the terms with subscripts are defined in Table I and Section 3.1. J is 
the unit matrix, H is a square matrix each element of which is unity, and p, £, 
are defined in (28). 


Solution of (31) yields the natural frequencies and associated elastic displace- 
ment modes. The absolute displacement modes are then obtained from (17), 
Section 6, with the use of (26)— (28), Section 6.1. 


8. Solution of Equations 


Once set up in numerical form the equations may be solved by any convenient 
method. It is suggested that matrix iteration’ or the escalator’? method is suitable. 


APPENDIX 
Treatment of Swept and Cranked Wings 


It has been assumed in the analysis that the wing and tail have straight 
flexural axes of the conventional type, but in modern aircraft this concept of elastic 
properties is frequently not even a good approximation to the true state of affairs. 
In particular, swept and cranked wings do not satisfy the conditions implicit in the 
assumption of a straight flexural axis. If it is assumed that chordwise sections of the 
wing do not deform (an assumption which is made for stressing wings with chordwise 
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Fig. 5. 


Swept wing representation. 


ribs), the problems which arise are related only to the calculation of influence 
coefficients, while the essential form of the dynamical problem is in no way altered. 
In fact, if influence coefficients are measured by static test of the structure, flexural 
axis considerations are irrelevant. 


For a wing of this type the idealised mass distribution is selected in exactly the 
same way as before, but no simple physical form is assumed for the elastic structure. 
A swept wing representation may appear as in Fig. 5. When a flexural axis exists 
it is assumed that the wing root forces and moments are transmitted to the fuselage 
at the point of intersection of the flexural axis and the fuselage centre-line. This is 
largely a matter of convenience since the idealised point of attachment of wing and 
fuselage is to some extent arbitrary; and similarly in the more general case the wing- 
fuselage junction point may be selected at some point within the root chord. This 
point then becomes the reference station on the fuselage. 


It now remains to make the appropriate modifications to D of (14) to cover 
the more general elastic properties. To do this it is necessary merely to define 
certain components of sub-matrices in terms of their physical significance and to 
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TABLE II 
Location in D | oor 
Modification 
Row | Column 
4 | 4 | Rw) 
| | _is replaced by matrix whose 
| | elements a, denote linear displacement of i“ mass of 
| wing due to unit linear force at j‘ mass (root encastré). 
4 | 7 (- R, 
is replaced by matrix whose 
elements a,,; denote linear displacement of i mass of 
wing due to unit torque at j“ spanwise station (root 
encastré). 
7 | 4 Rw) 
| is replaced by matrix whose 
| elements a; denote angular displacement at i span- | 
| wise wing station due to unit linear force at j” wing | 
| mass (root encastré) (transposed of above). 
7 | 7 


is replaced by matrix whose | 
elements a, denote angular displacement at i” wing 
station due to unit torque at j™ station (root encastré). 


replace them by the corresponding matrices thus defined. The appropriate 
modifications are set out in Table II*. 


Similar modifications may be necessary for matrices relating to the tail. The 
above table applies similarly to the tail when the subscript “t” replaces “w,” 
“tail” replaces “ wing” and the row and column numbers are ordered as follows: 
2,2; 2,6; 6,2; 6,6. 
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Airscrews' at Supersonic Forward Speeds 


J. C. BURNS 
(Department of Mathematics, The University, Manchester) 


(Summary: The flow past an airscrew rotating with uniform angular velocity in a 
uniform supersonic stream is considered from the point of view of axes fixed in the 
airscrew, the motion relative to these axes being steady. A linearised equation and 
boundary conditions are found for the potential of the disturbance flow caused by 
airscrews which produce only small perturbations in the main stream. This equation 
is solved by expanding the potential in powers of the ratio, ,, of the tip speed of the 
airscrew to the speed of sound in the undisturbed stream. Equations and boundary 
conditions are found for the coefficients of each term of the series and it is seen 
that the term independent of x satisfies the ordinary linearised potential equation for 
fixed axes. By subtracting suitable special integrals, the more complicated equations 
for the coefficients of the various powers of x can be reduced to this ordinary potential 
equation. Hadamard’s methods have been applied by Puckett’, Evvard® and, 
finally, Ward" to the problem of finding the flow past a thin wing fixed in a 
supersonic stream and their methods are applied to find the terms of the series in 
our case. The first two terms are found explicitly for a particular type of windmill. 
Approximate expressions are then developed for the torque and the drag on such 
a windmill and the results applied in special cases. The best working conditions 
for windmills of this type are found by considering the effect on efficiency and 
power of varying the shape and the speed of rotation, within the limits imposed by 
the strength of the material of the airscrew.) 


1. Introduction 


In the supersonic motion of an aircraft or projectile through the air, it may be 
desirable, in some cases, to have a source of power other than the main driving 
force (if any). An airscrew mounted on the body and designed so that the 
aerodynamic reaction of the air produces a torque which tends to increase its 
angular velocity (i.e. a windmill) could be used as such a subsidiary power source. 
The present paper describes the flow past airscrews of this type. The theory which 
is developed could be modified, without any essential change, to deal with propellers; 
the specialisation to the case of windmills rather than to propellers is suggested 


+The word airscrew is used in this paper in its general sense, referring to any type of screw 
designed to rotate in air. Although the theory can be modified to deal with propellers, 
it is applied here to windmills.—Eb. 
Paper received November 1950. 

[The Aeronautical Quarterly, Vol. III, May 1951] 
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by the greater possibility of practical application and does not relate to any 
limitation on the theory. 

Power is obtained from the windmill at the expense of the forward momentum 
of the body, the reaction of the air producing a drag as well as a torque on the 
airscrew, and the efficiency of the windmill is measured as the ratio of the power 
produced to that required to overcome the drag. For high efficiency, it is desirable 
that the drag be as small as possible. The airscrews considered are thus taken to be 
thin aerofoils approximating to that aerofoil of zero thickness which, when rotating 
at the operational speed of the airscrew, produces no disturbance at all and therefore 
suffers no drag (or torque). Airscrews of this type produce only small disturbances 
in the fluid and it is assumed that the squares and higher powers of the perturbation 
velocities may be neglected. This means that linearised equations of motion may 
be developed and that entropy changes across any shock waves in the fluid may be 
neglected. It is assumed also that the effects of viscosity in the flow are confined 
to a negligibly thin boundary layer on the surface of the airscrew so that the flow 
may be considered to be isentropic. Finally, it is assumed that the air may be 
treated as a perfect gas, which is an excellent approximation under a wide range of 
conditions, including those in the atmosphere. 

Airscrew theory in the subsonic case has been studied in great detail and an 
account of the theories which have been developed is given, for example, in Ref. 4 
from which the following brief outline is, in the main, derived. References to the 
authors mentioned here, and many others, are given by Durand. The earliest 
theory, the “ momentum theory” of Rankine and Froude, ignored the geometrical 
shape of the airscrew, replacing it by an “actuator disc” and deriving the resultant 
forces on the airscrew by consideration of the overall changes in the momentum 
of the fluid. The “blade element” theory, due initially to Froude and mainly to 
Drzewiecki and Lanchester did treat in detail the flow past the airscrew by dividing 
the blade into elements along its radius and treating each of these as an aerofoil 
section moving in a main stream which is modified to allow for the transverse 
component of velocity caused by the rotation. It was found that this theory could 
explain the behaviour of an airscrew in general terms only and that it failed to 
give accurate numerical results. Moreover, the theory suggested that the efficiency 
of an airscrew tends to unity as the drag tends to zero, a serious divergence from 
the results of the momentum theory which set an upper limit to the efficiency. 


Side by side with the blade element theory for airscrews, there has developed 
the “ quasi-static” theory for treating the flow past wings in roll. In this theory, 
allowance is made for the (small) rotation of the wing, not by modifying the main 
stream as in the blade element theory but by assuming a suitable change in the 
shape of the wing section (more accurately, by making an appropriate change to 
the boundary condition at the surface of the aerofoil). It is apparent that these two 
theories achieve the same results for the wing in roll. Although the quasi-static 
theory has the advantage in this case that it enables the problem to be solved by using 
the ordinary methods for finding the flow past an aerofoil at rest in a uniform stream, 
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the blade element theory is more general, being applicable without the restriction to 
slow rates of rotation. 


Eventually, the key to the completion (and reconciliation) of the theories was 
found by Lanchester to lie in the vortices shed by the airscrew and composing its 
wake. These produce an interference velocity at the blades which has to be taken 
into consideration when modifying the incident velocity at the surface of the aerofoil. 
The difficulty of the “ vortex theory” is in the analysis of the vortex distribution in 
the wake and the derivation of the correct contribution to the incident flow. 
Betz, Prandtl and Goldstein attacked the problem by considering the actual helical 
vortex sheets in the wake but a common procedure has been to simplify the problem 
by assuming that there is a large number of blades and postulating a uniform 
distribution of vorticity in the wake, thus neglecting the periodicity associated 
with a small number of blades. 

In supersonic flow, the problem is much easier since the effects of disturbances 
in the fluid are confined to a region downstream from the points where they occur, 
so that the trailing vortices have no influence whatever on conditions on the surface 
of the blades of the airscrew. For this reason, too, there is no interference between 
the blades of the airscrew, provided that each blade lies outside the region of 
influence of the others. For high pitch, which would be necessary at high speeds 
to keep the angular velocity down to a level which the materials can stand, this 
condition is clearly satisfied. In the supersonic case then, the method of the 
blade element theory, or (as is preferable) of the quasi-static theory, may be used. 
However, it is now possible to go further and demonstrate that the solution obtained 
from this theory is the first approximation for a slowly rotating airscrew and to find 


the second and higher approximations. 

The problem is considered from the point of view of axes fixed in the airscrew 
which is taken, for convenience, to be rotating with uniform angular velocity 2 
about a fixed point in a uniform supersonic stream of velocity U parallel to the 
axis of rotation. The linearised equation for the disturbance potential is found and 
solved approximately for airscrews for which the ratiot x of the tip speed 2b to 
the speed of sound in the main stream is small. The potential is expanded as a 
power series in x and a method given for finding the terms of the series in succession. 
In each case, the problem is reduced to that of finding the flow past an aerofoil fixed 
in a uniform supersonic stream using the methods developed by Puckett"), Evvard"’ 
and Ward". The first term of the series is exactly the solution derived from the 
quasi-static theory and the later terms give the higher approximations to the flow. 

The theory is applied to a particular type of airscrew for which the first two 
terms of the series for the potential are worked out. Approximate expressions are 
then found for the drag and the torque on an airscrew of this type and the results 


ty=2b/a,=MX=7M/J where M=U/a,. X=b/U is the parameter defined by Durand 
(Ref. 4, p. 328) for windmills and J is the advance/diameter ratio defined for propellers by 
Durand (Ref. 4, p. 173). 
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applied in special cases. The power which can be derived from such airscrews and 
the efficiency at which they work are found for a range of values of a parameter 
which determines the shape of an airscrew of given plan form and thickness, 
rotating at given angular velocity. It is found that the efficiency reaches a maximum 
for a certain value of this parameter which therefore gives the best working 
conditions. The variation in this maximum efficiency and the corresponding power 
with the angular velocity is then considered. Best results are obtained for high 
values of the angular velocity, provided that the tip speed does not exceed a limit 
imposed by several restrictions, some of a physical nature and others concerned 
with the limitations of the linearised theory. Finally, an estimate is made of the 
way in which the strength of the material of the airscrew restricts the shape. 


Notation 


a, B=x +By } characteristic co-ordinates for the equations for 
a’, B’=x'+ By’ potential flow relative to fixed axes 


¢ parameter determining the thickness of the airscrew 


(=P/(xc*p,U*) power coefficient 


< 


(m Value of ¢ when airscrew is working at maximum 
efficiency 

k,¥ coefficients depending on mw and x and defined by 

n=(A+K 


A 

Am Value of A which gives maximum efficiency 
nh bie 
Seé Kk, Vv 


p density of air 
p, density of air in undisturbed main stream 
@ perturbation potential 
%(:)s9%(2) two parts of as defined in Section 3 
Pos Pas - } 
Poti» Pili)» + 


$,,%.,. - solutions of the equations for potential flow relative to 
fixed axes associated with 9,, ¢., . 


x=Qb/a, ratio of tip speed of airscrew to speed of sound in main 
stream 


© angular velocity of airscrew 


coefficients of powers of x in series expansions of 9, $,) 


defined by the relation 2. + =angular velocity at which 
mean surface of airscrew causes zero disturbance 
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ind a 
ter 
SS, 
1m b 
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yer 
gh 
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h 
F 
I (2, B) 
J 
or 
M 
n 
n n* 
P 
y 
p 
Ap 
Po 
rr 
S 
U 
Vv 
Ww 
X 
4 \ 
x’, y’, 2’ 


local speed of sound 

speed of sound in undisturbed main stream 

radius of airscrew 

(M? — 1)! 

chord of airscrew 

drag 

Young’s modulus 

aerodynamic force on airscrew 

quadratic moment of root section of airscrew blade 
integral in Section 3.11 


advance/diameter ratio for propeller—see footnote, 
p. 25 


Mach number of undisturbed main stream (=U /a,) 
normal to surface z=y tan (0 x/U) 

normal to actual airscrew surface 

power derived from airscrew (=2 Q) 

airscrew torque 

pressure 

perturbation pressure 

pressure in undisturbed main stream 

position vectors 

surface of integration in Section 4 

maximum stress in airscrew blade in Section 4.1 
velocity of main stream relative to fixed axes 
velocity of main stream relative to moving axes 
perturbation velocity 

stream velocity relative to fixed axes 

stream velocity relative to moving axes 

perturbation mass flux vector 

ratio of tip speed of airscrew to main stream velocity 


(=2b/U), see footnote, p. 25 


components of position vectors 
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2. Mathematical Formulation 


Let x, y, z be a right handed system of Cartesian co-ordinates such that x is 
measured in the direction of the main stream, which has velocity U relative to the 
axes. Let the axes be fixed in the airscrew so that they are rotating with constant 
angular velocity €2=(— ©, 0, 0), the angular velocity being positive in the direction 
of forward motion of the airscrew. Since the airscrew is constantly at rest relative 
to these rotating axes, the motion referred to them is steady. 


Let v be the fluid velocity relative to fixed axes with which the rotating axes 
instantaneously coincide and let v’ be the velocity relative to the moving axes; then 


If the motion is assumed to be irrotational relative to the fixed axes, then 
curlv=0._ . : (2) 


Bernoulli’s equation for the pressure in steady motion is (Ref. 5, p. 20) 


v?—v. Qart+ | dp/p=constant (3) 


Nw 


The equation of continuity in steady motion is (Ref. 5, p. 12) 
div (pv’)=0. é (4) 


Equations (1) to (4) are the equations of motion for the steady irrotational flow 
past the airscrew. 


Let v=U+u . . (5) 
where, by (2), u=grad¢; ‘ . (6) 


@ is the perturbation velocity potential. It is assumed that the disturbances caused 
by the airscrew are sufficiently small for terms quadratic in ¢ to be neglected. 


Let p,, p, and a, be the pressure, density and speed of sound in the undisturbed 
stream (where v=U). Then, from Bernoulli’s equation (3), 


p 2 
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is the main stream velocity relative to the rotating axes. Equations (1), (5) and 
(8) show that 


By expanding the integral on the left hand side of (7), the pressure and density are 
obtained in terms of the velocity in the forms 


and p=p,(1-U’.u/a,?)+O \ul?. 


The vector w can be interpreted as a perturbation mass flux vector in the sense that 


Finally, substituting for pv from equation (13) in the continuity equation (4) gives, 
correct to O( |u| ), 


By use of equations (6), (8) and (11), equation (14) becomes the equation for the 
disturbance potential , namely 


2 2 2 2 2 2 


0z0x Ox ay 
(15) 
where M=U/a,. 
In cylindrical polar co-ordinates x, r, 9 for which 
y=rcos@, 
the equation is 
109 1079 07079 20U ao 
(l M53 + oF ror a,? @xd6 (16) 


Equations (15) and (16) will be satisfied by the disturbance potential of the flow past 
an airscrew rotating with angular velocity (—, 0, 0) in a uniform main stream with 
velocity (U, 0, 0), provided the perturbations are small. It has already been pointed 
out that for this to be so, the airscrew must approximate to the surface (of zero 
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thickness) which would produce no disturbance at all when rotating in the same 
way. The normal nto such a surface must satisfy the condition 


where U’ is given by (8) as the velocity of the undisturbed main stream relative to 
the axes fixed in the rotating surface. If the surface has equation f(x, y,z)=0 
relative to these axes, condition (17) gives as the equation for f (x, y, z): 


The surface which satisfies this equation and passes through the line x=z=0 
has equation 


and the equation of the airscrew surface is assumed to approximate to (19). 


The boundary condition on the surface of the airscrew is, exactly as in (17), 
that the normal component of the fluid velocity relative to the airscrew be zero, 


where m* is the normal to the surface and, since the surface differs only slightly 
from (19), can be taken as differing from m by a vector whose magnitude is O |ul. 
Substituting for vw from equation (9) gives finally as the fully linearised boundary 
condition that 


which, to a first order approximation, may be applied on the surface z= y tan (Qx/U) 
over the plan form of the airscrew. 


The complete mathematical problem to be solved is given by equations (15) 
or (16) under boundary condition (21), subject to the vanishing of the disturbance 
on any plane, normal to the flow, which lies ahead of the airscrew. 


3. Solution of the Potential Equation in Series 


The equations and boundary conditions formulated in Section 2 apply to all 
airscrews, whether propellers or windmills, but in developing the method of solution, 
it is convenient to consider a particular type of airscrew and, for the reasons 
indicated in the introduction, a windmill is chosen. 


If an airscrew of zero thickness whose surface is given by equation (19) rotates 
at angular velocity ©, it causes no disturbance in the air and there are no forces 
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acting on it. If now the speed of rotation is reduced slightly, there will be small 
disturbances in the air and the reaction of the air will produce a torque which will 
tend to restore the angular velocity to the value 2. If this torque is utilised in doing 
work and the velocity maintained at a value less than ©, the airscrew is being used 
as a windmill to provide power, the amount of which depends on the extent to which 
the angular velocity is reduced. 


Accordingly, an airscrew is considered which has its mean surface given by 


(where » is small), but which rotates at speed ©, so that there will be a torque 
acting on it and tending to increase the angular velocity. It is necessary to give 
the airscrew a finite thickness and the upper and lower surfaces of the airscrew 
are taken to be given by the equations, 


(where < is small). For convenience, the symmetrical case only is considered but 
it may be seen that no difficulty of principle arises in the unsymmetrical case. 
If propellers rather than windmills were being discussed, the equations of the 
surfaces of the propeller would be the same as (23), except that (0+) would be 
replaced by (Q —w). 


The potential equation (15) is now solved, with boundary condition (21), for 
the airscrew given by equation (23). Consider first the upper surface (23) and 
suppose that squares and higher powers of <« and » may be neglected. Applying 
boundary condition (21) to this surface will lead to a condition of the form 


=» tan © g, (xX, y) +e g, (x,y). (24) 

Similarly, applying (21) to the lower surface (23) leads to the condition 
(0 /An)z=y tan x/U)-0 = OB, (xX, y)— (x, y) (25) 
@ can thus be split into two parts ¢,,) and 9.) which satisfy boundary conditions: 
(0,1) / On)z=y tan tan Ox/U)-0=O8, (x,y) (26) 
and (06,2) / On)z=y tan x/U)+0= — (09()/ON)z=y tan Ox/U)-0= (27) 


It is easily seen that ¢,,) is the potential for the flow past the aerofoil whose equation 
is obtained by putting « zero in equation (23) and is thus given by (22). Similarly, 
%(2) is the potential for the flow past the aerofoil whose equation is obtained by 
putting » zero in (23) and so is 


Methods for finding ¢,,, and ¢,,) are given in Sections 3.1 and 3.2. 
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3.1. Atrscrew of zero thickness designed to rotate at angular 
velocity (+) without disturbing the flow 


In finding the first part %,,) of the potential for the flow past the airscrew with 
surfaces given by equations (23), the flow must be found past the aerofoil of zero 
thickness which rotates at angular velocity 2 but is designed to rotate without 
disturbing the flow at angular velocity (2+) and so is given by equation (22). 


In theory, this problem can be solved directly by means of Hadamard’s 
general theory for boundary value problems for hyperbolic equations, using methods 
analogous to those for the ordinary problem of an aerofoil fixed in a uniform stream. 
Attempts to do this were unsuccessful because of the complexity of the integrals 
involved. However if the potential ¢ is expanded in powers of the ratio of the tip 
speed, 2 b, to the speed of sound in the main stream, a,, as 


oot ...., (29) 


Hadamard’s method can be applied to find each of the terms 9,, 9,,. . . . in turn. 
In fact, the problem to be solved in each case reduces to solving the equation 


(1 — M?) 079 / + dy? + ‘ 


for the flow past an aerofoil fixed in a uniform main stream. Hadamard’s method 
has been used to solve this problem in detail and a complete account of this is 
given by Ward". 


The procedure is to substitute for from (29) in the equation (15) and boundary 
condition (21) and, by equating coefficients of powers of 2b/a,, to obtain the 
equations and boundary conditions satisfied by 9, %,,... This leads to the 
following equations for the terms of the series 


(33) 


Only 9%, and ¢, are considered in what follows, and, for determining the 
boundary conditions on these, and » may be assumed sufficiently small for w? 
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and 9° (and hence, also, for example, 2? and higher powers) to be neglected. 
In this case, the direction cosines of the normal m* to the surface (22) are 


~(Q+o)y/U, —(Q+o)x/U, (x? +y?)/(2U2). (34) 

It follows that the boundary condition (21) becomes 
0¢/On=woy onz=ytan(Qx/U) . (35) 
to the order of the approximation and, as (26) indicates, this condition holds on 


both sides of the surface. An expansion of (09/0), _, means of Taylor's 
theorem gives 


a a Q 

Co Co 
(22) 
z=y tan (Ox/U) On/ U CZON’ z=0 


to the order of the approximation previously defined; and since 09/dn is the 
component of grad ¢ along the normal to the surface z=y tan(x/U), 


Substituting for @ from (29) and combining the results of (36) and (37), from the 
boundary condition (35) are obtained the following boundary conditions on ¢, and 
¢,, over the plan form: 


261) 0o ) 
To obtain the conditions on ¢,, ?,, . . . it would be necessary to retain terms in the 


cube and higher powers of 2 in forming condition (35). The condition on ¢,, for 
example, is 


7 
2xy dz Oydz 2 d2° 
(40) 


xy 


The equation (31) and boundary condition (38) satisfied by 9%, are identical 
with those given for the potential by the quasi-static theory. On this theory the 
airscrew is treated as an aerofoil fixed in a steady stream and the (small) effect of 
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the rotation is allowed for by amending the boundary condition at the surface of the 
airscrew. This means that the potential is given by equation (30), the ordinary 
linearised equation in fixed axes, but the boundary condition becomes 


(09/0z).-0o=(—U’ . m*):=0; . (41 


in the case of a fixed aerofoil in a uniform stream, the right hand side of (41) would 
be (—U . m*)._». For the particular airscrew under consideration, (41) reduces to 


(0¢/0z):-o=wy, 


which is the same as condition (38) for ¢,. Thus the quasi-static theory gives the 
first approximation to the flow past the airscrew for small angular velocities. The 
method has been used, for example by Robinson’, in the case of a delta wing. 

The form of the boundary condition (38) on 0¢,/0z shows that 9, is an odd 
function of z and, since equation (31) is the usual linearised potential equation, 
@, can be found by the method given by Ward? for the anti-symmetrical part of the 
potential for the flow past a thin aerofoil fixed in a uniform supersonic stream. 

The first step in solving equation (32) for ¢, with boundary condition (39) is to 
find a particular integral of the equation. Using this integral it is then possible to 
reduce the problem again to the solution of the ordinary equation for a fixed 
aerofoil, this time in the symmetrical case treated by Puckett" (and also given by 
Ward’). In cylindrical polar co-ordinates, equation (32) becomes 


and the expression 
Mx 0, Mx 0%, ty) 


can be shown, by substituting in the equation, to be a particular integral. Thus 


__ Mx z 


9,=9,+ 


provided that 4, satisfies the equation 


(1 —M?’) 


Substituting for ¢, from (44) in boundary condition (39) gives as the condition on 
é, over the plan form: 
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Since ¢, is an odd function of z, this condition is such that 
/0z). +0- (0¢, /0z): 0 (47) 


so that 4, is an even function of z and can be found by Puckett’s method. 


In the same way, the problem of finding 9,, 9,,. . . can be reduced, by finding 
particular integrals of the equations, to that of finding, alternately, the symmetrical 
or the anti-symmetrical part of the potential past a fixed aerofoil in a uniform 
supersonic stream. For example, 


Mx 06, Mx? 0° x 
bB? 30 +4 Bw) 597 . (48) 
0 
where B?=M?-1 
25, 06, 76, 
and 


and 06,/éz is given by substituting for 9, from (48) in condition (40) for 
(0 92/02), _ 


3.1.1. Solution for ¢, for an aerofoil of rectangular plan form and independent 
wing tips 


@, Satisfies equation (31): 


and boundary condition (38) on the plan form: 


Solving for ¢, in the particular case of rectangular plan form, let ABCD (Fig. 1) 
be the plan form with span AD=25 and chord AB=c. The origin of axes is the 
mid-point of the leading edge with the x-axis in the direction of the stream and 
the y-axis along the leading edge. Let AF and DE be the characteristic lines for 
equation (31) through A and D respectively and lying inside the plan form. 


Then, from Ward’s’ equation (13), at a point P(x, y,0) inside AFED, 


y dx’ dy’ 
0)= — | { — } 
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x 
Fig. 1. 
Plan form of airscrew. 


where the integration is taken over the area of the plane z=0 for which 
x =< x+B(y-y). Following Ward, the integral in (49) is transformed by 
introducing characteristic co-ordinates (a method due to Evvard): 


and the equation 


_ o (2’ — B’) da’ dB’ 

is obtained, the integration being taken over the area 2’< 2, #’< 8. The integral in 
(51) is split into two parts so that 


6) 
= [7 (z, B) I (8, 2) | (52) 
4 
where I (2, 8)= | 
Then, from (52) and (53), 


(54) 


A y 
/ 
/ F 

\ / 
\ / 
\ / 
\ / 
\ / 
\ / 
\ 
\ / . 
: 
\ 
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Again, for a point P(x, y,0) on the wing tip ABF, from Ward’s® equation (30), 
», is given by an integral of the same form as (51) but over the area KLPM in 
Fig. 2. The integral can be split into three parts and eventually gives 


(x, y,0= = { Bxysin x 
x [B(b—y)(x+ B(y—b))} } 
(55) 


A similar result can be obtained for the wing tip CDE. 


@, has thus been found at all points on the surface of the aerofoil. It may be 
noted that 9, and @9,/0@x are both continuous across the line AF and across the 
wing tip AB (9, =0 for |y|>5). 


3.1.2. Solution for ¢, for an aerofoil with rectangular plan form and independent 
wing tips 


As has already been shown, the problem of finding ¢, reduces to that of finding 
6, which is given by equation (45): 


with boundary condition (46) on the plan form: 


¢ ~ Mb oz 
Since 6, is an even function of z, 
(04, /0z):-0=0 


at points of the plane z=0 outside the plan form, and 4, (x, y,0) is given, at all 
points of the plane, by 


(56) 


the integration being over the area for which 


x = x+Biy-y’. 
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(«,-«) [(e-288),- (2-288) 
/ 
a= 
/ — 
/ 
/ 
/ 
B 
Fig. 2. 
Now, for points in AFED, 
(06, /0z).-0= — (x? — B* y*)/(Mb (57) 
so that for points P(x, y,0) in AFED, 
or, in characteristic co-ordinates, 
= 48 Mb —1448+5 B?). (58) 
Hence (x, y> 0)= - 6 B’ y’ ) (59) 
and 9, (x, y,0)= (x? + 6 Bt (60) 
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In the same way, using the value of , given by (55) it is possible to find an 
expression for ¢, and hence for ¢,, for points in AFB. The integration is not 
completed here but it is easily seen that the value of ¢, obtained will be continuous 
across the characteristic AF with the value given in (60). However, for points 
on the edge AB, where there will in fact be high velocities, the linearised theory 
breaks down and ¢, becomes infinite owing to the presence in 0 4,/0z (and 2 4,/0z) 
of the term in 0? ¢,/0z?. 


3.2. Symmetnical airscrew 


The second part, ¢,,.), of the potential for the flow past the airscrew given by 
equation (23) is found by considering the flow past the aerofoil given by 
equation (28): 


z=ytan(Qx/U)+-f (x,y). 
If £2 and < are assumed small so that <* and {2° (and hence also (2? < and higher 


powers) may be neglected, the direction cosines of the normal m* to the upper 
surface given by equation (28) are 


Boundary condition (21) then becomes 


(0¢/ On); y tan (Qx/U)+0 = Us Of /Ox ° ° (62) 


and, in the same way, considering the lower surface given by (28), the boundary 
condition 


(0¢/ On). y tan (Qx/U)-0= — Use of /Ox (63) 


is obtained. Expanding ¢ as the power series given by (29) and substituting in 
(62) and (63) gives as the condition on 9, 


/0z)2=- +0= . . (64) 
while the conditions on ¢,, ¢, . . . and the equations for ¢,, %,, %,.. . are exactly 


as given in Section 3.1, equations (39), (40) and (31)-(33). 


In this case, the boundary condition (64) on ¢ ¢,/0z is such that 9%, is an even 
function of z and so can be found by the use of the method given by Puckett" for 
the symmetrical part of the flow past an aerofoil in a uniform stream. 9, and the 
associated 4, will now be odd functions of z so that the method required is that 
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given by Ward"? for the anti-symmetrical part of the flow past an aerofoil in a 
uniform stream. 


3.2.1. Solution for ¢, and , for an airscrew with symmetrical circular arc cross 
section 


As an example of the application of the method, an airscrew is considered 
which is given by equation (23) with f (x, y)=4 (cx —x*)/c. This gives a cross section 
of two equal circular arcs with end points x=0 and x=c and maximum 
thickness 2 < c. 


The potential 9, for points independent of the wing tips is simply the well- 
known potential for steady two-dimensional flow past a circular arc aerofoil and 
is given by 


{ (x— Bz}? —c(x—Bz)}. (65) 


With this value for ¢,, the boundary condition on the function ¢, associated with 
the second approximation 9, becomes 


(0¢,/0z):-0=[4U y/(M Bbc)] (4x—-c). . (66) 
Hence, using methods similar to those in Section 3, it is found that 

@, (x, y,0)=[4 U xy/(M B? bc)] (c—2 x), 
and hence that 

(x, y,0)=[4 U < xy/(M bc)] (2 . (68) 


for points not affected by the wing tips. 


3.3. Summary of results 


By superposing the results of Section 3.1 and 3.2, the following solutions are 
obtained, at points on the part of the surface of the airscrew independent of the 
tips, for the first and second terms of the series for the potential for the flow past 
the airscrew surface given by the equation 


z=ytan[(0+0)x/U]+4¢ (cx—x*)/c ; (69) 

namely, (x, y,0)= —o xy/B+4U —cx)/(Be) . (70) 
Wx 

and 1 (%¥,0)= (x? + 6 BY y*) + Mie 
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The values of 9, and ¢, have been found at points on the surface of the airscrew 
because these are the values required in calculating the forces on the airscrew. 
The methods of Ward and Puckett do however give the values at points off the 
surface of the airscrew without difficulty. For example, the expression for ¢, for 
points which are not affected by the wing tips is (for z > 0) 


(x, y, — wy (x/B-z)+ [(x — —c (x— Bz)]. 


(On putting z=0, this result reduces to equation (70) ). 


When the results of this section are used in Section 4 for calculating the 
aerodynamic forces on an airscrew the comparative magnitudes of the two small 
quantities « and , which have hitherto been quite independent of each other, will 
have to be taken into consideration. In fact, it will be noticed immediately that < 
always appears in these results multiplied by the large velocity U so that it is to 
be expected, for example that terms of the first order in < will be larger than those 
of the first order in w. 


4. Aerodynamic Forces on the Airscrew 


The principle of the conservation of momentum can be expressed, in moving 
axes, in the form 


{ Apn+pv(v'.m)} dS=0 ; 
where S is a closed surface fixed relative to the moving axes, \p is the perturbation 


pressure and v, v are the fluid velocities relative to the fixed and moving axes 
respectively. 


Thus, if S is a closed surface surrounding the aerofoil and sufficiently far from 
all singularities for the flow to be given accurately by the linearised potential and 
for |V¢| to be uniformly small on it, the aerodynamic force F on the airscrew will be 


F=— | { dS. 
AY 


With the aid of the continuity equation (4) and the fact that pw’=p,(U’+w) 
(equation (13)), equation (74) becomes 


F=- | { Apn+p, grad ¢(U’+w).m} dS >) 
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where Ap and w are given by equations (10) and (11). By an argument exactly 
similar to that used by Lighthill' for the case of fixed axes, it can be shown that 
when F is given by the integral in (75), the surface S can be taken to be any closed 
surface surrounding the airscrew. 


Suppose the airscrew given by equation (23) has a plan form symmetrical about 
the axis of rotation, and that f(x,y) is an even function of y, as in the particular 
example given by equation (69). It has been shown in Section 3 that the potential 
for the flow past such an airscrew can be split into two parts, ¢,,) and %,,), and each 
of these has been expressed as a power series in 2.b/a,. In finding the forces on 
the airscrew, the symmetry properties of the various terms are important. It can 
easily be deduced from the form of the equations and the boundary conditions that 
and po (2), P2i2). - - - are even functions of z. Also, for the particular type of plan 
form and cross section chosen, - ANd ate odd functions 
of y and (1), ANd are even functions of y. 


If now (75) is integrated over the upper and lower surfaces of the airscrew and 
use is made of these symmetry properties, the following expressions are obtained 
for the torque and the drag 


( a, U 


+O 1)? w 
(76) 


in the direction of the angular velocity and 


| ay 
p= 1) Poz (1) + Porte Poze y+ Qy + M? pe + 


+ Pi x(2) + toxin } aS + O(p,eos 2 c). 
(77) 


In each case the integration is over the half of the airscrew plan form for which y>0. 


Equations (76) and (77) show that the resultst obtained from the quasi-static 
theory for the torque and the drag, namely, 


and D= —4p, | | Poz(1) +0 dS (79) 


+Compare the result given by Ward®) (equation (60)) for the drag on an aerofoil. 
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are the first approximations obtained by retaining only the terms independent of ©. 
However, it will be seen when numerical values are inserted that, depending on the 
relative magnitudes of the physical quantities, these terms are not necessarily the 
largest ones. 


4.1. Application to special type of airscrew 

These results for the torque and the drag are now applied to the particular 
airscrew with surfaces given by equation (69) and with plan form which is rectangular 
but which has its corners cut off along the characteristics through the tips of the 
leading edge, i.e. the plan form is AFED in Fig. 1. By substituting in equations 
(76) and (77) the values obtained for %,;;), (2), %::2) from equations (70)-(72), 
the following results are derived for the torque and the drag. 


4 


(80) 
and 


32p,U2e 


3B 
(81) 
Consider the special case in which the main stream has Mach number M= 2 
so that B=1 and the velocity U=a,/2 where, as an example, the speed of sound, 
a,, and the density p, are assumed to have the values given by Pankhurst and Conn“? 
for conditions in the standard atmosphere at a height of 20,000 ft. 


If A, u, x and ¢ are non-dimensional parameters defined by the relations 


then the expressions for the power produced by the windmill and the power 
absorbed by the drag become 


Q=(V72/3)p.a, c* x x { (u—1)*] A+84)- 


3A 
[ue (u (u+4)] T o( \ wats , 
(83) 
and UD=(¥ 2/3) py 2? c* x 10-* x { 
+ “(2u-1)4+0G: watts . 
(84) 
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From (83) and (84) it can be seen that the efficiency, »=(2Q/(U D), depends 
only on the four parameters A, », ¢, x. For given values of » and x, the efficiency 
takes the form 


A+Ké 
= 


where « and v depend only on « and y. The way in which » varies with A and 
depends on the magnitudes of « and v and it may be seen that there are three main 
types of behaviour possible. However, with the range of values of « and v required 
by the approximations of the theory, only two of these types occur. In each of these, 
n is bounded for positive values of A and reaches a maximum value, »,,, for a certain 
value, A», of A, with corresponding value, Gn, of the “ power coefficient” (. A, is 
given by the equation 


Aw { (k?+v)e?—Ke }4 KE, (86) 


It is clearly desirable to operate the windmill at the highest possible efficiency; 
therefore, keeping (2c fixed, the variation with » and « of A,, is considered and the 
corresponding values of »,, and (,. This variation is shown in Figs. 3, 4 and 5, 
for values of » between 2.0 and 10.0 and of « between 0.01 and 0.05, keeping 
Qc=1,000 cm./sec. Fig. 5 also indicates the actual power obtained from a 
windmill of chord 2 cm. under these conditions. From the results shown in these 
three diagrams, it is seen that the best working conditions for given 2c are obtained 
by taking < as small and » as large as possible. This means taking the thickness 
of the airscrew as small and the span as large as possible for given chord and, 
in each case, the limit will be imposed by the strength of the material of the airscrew. 
These conditions not only give the highest value for »,, but also ensure that A,, has 
the lowest value possible. This is desirable as the linearised theory has been 
developed on the assumption that » and hence A is small. 


To obtain an indication of the limits to u and ¢, an estimate is made of the 
maximum stress in the blade and this is compared with the elastic limit for the 
material of the airscrew. The bending moment due to the pressure forces on the 
blade is a maximum at the root of the blade and is equal there to half the total 
torque on the airscrew. The greatest thickness of the blade at the root is 2¢c and 
hence, by the well-known formula relating the stress in a beam to the bending 
moment (see, for example, Lamb"’), the maximum stress, S’, is given roughly by 


where 4 Q=bending moment at root of blade 


and I=quadratic moment of the root section with respect to its neutral line. 
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Fig. 3. 
Variation of A,, with » and ¢ for 2c = 1,000 cm. /sec. 


The maximum stress must be less than, say, half the elastic limit of the material 
of the airscrew. In Fig. 6 are shown the values of the ratio of these two stresses 
for different values of » and < when the airscrew is working at maximum efficiency 
and c= 1,000 cm./sec. It is assumed that the airscrew is made of spring steel with 
an elastic limit of 70,000 Ib./in.2 From these results, the values of » and « can be 
obtained for which the maximum stress is less than half the elastic limit (taking 
a safety factor of 2) and it is seen, for example, that « must not be greater than 3.0 
when ¢=0.03. In practice it would be necessary to consider sets of results of this 
type for different values of the angular velocity. 


The variations in A,,, 7 and ¢,, with y, keeping « and « constant are now 
considered. In Fig. 7, the values of these parameters are given for values of x 
between 0 and 0.76, keeping «= 3 and «=0.03, as a particular example. For values 
of x greater than approximately 0.93 the variation of » with A changes over to the 
type in which » decreases steadily from infinity to zero as A goes from zero to 
infinity. This is a non-physical solution and is ruled out because these values of 
x are evidently too large for the linearised theory to be valid. The variation in ¢. 
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Variation of »,, with » and ¢ for Qc= 1,000 cm./sec. 


with x also shows how the magnitude of x must be limited. The initial steady 
increase of (¢,, with x is followed by a rapid decrease after the maximum is reached. 
It can be seen from the algebraic dependence of A,, and ¢,, on x that they each become 
zero for values of x approximately equal to 0.93 while, at the same value of x, 
4 is approximately unity. This decrease in power with increase of tip speed is not 
to be expected on physical grounds and indicates that the linearised theory is not 
valid for such large values of x. That this is so is further indicated by the fact that 
there is no significant change in the shape of the airscrew at or near the value 
x=0.93. The equation (19) of the basic airscrew surface can be written 


z=ytan[y x/(u Mc)] . (88) 
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Variation of ¢,, and power 2Q (for c=2 cm.) with « and ¢ for 2c=1,000 cm./sec. 
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Variation with » and e of the ratio of the maximum stress in the 
airscrew blade to the elastic limit of steel. 
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Variation of A,, and with x for «=3, «=0.03. 


Tm 
which, for »=3, M=+/2 and x=0.93, becomes approximately z= y tan (13 x°/c). 
It is evident that this surface is twisted only slightly from the plane z=0 and does 
not change in any important way when x changes slightly. In the case considered, 
it seems reasonable to suppose that the theory remains valid for values of , up to 
0.5. As far as variations of A,, and 7, with , are concerned, the best results are 
obtained for high values of x which produce high values of »,, and low values of A., 
the merits of which have been discussed previously. 


5. Conclusion 


Apart from the estimate of the limits imposed by the strength of the material 
of the airscrew on the aspect and thickness ratios of the blades, no attempt has been 
made here to effect the necessary compromise between theoretical and practical 
requirements in designing a supersonic windmill. An indication of the procedure 
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Fig. 7(b). 


Variation of ¢,, and power 2Q (for c=2 cm.) with x for »«=3 and e=0-03. 


can be obtained from the results which have been given in Section 4. To take a 
specific case, it may be seen from Fig. 6 that the values »= 3, «=0.03 appear to give 
a practicable shape for a steel airscrew of the type considered, working at maximum 
efficiency and rotating so that Qc=1,000 cm./sec. (and ,=0.1 approximately). 
Reference to Fig. 4 shows that, under these conditions, this particular windmill has 
an efficiency of about 20 per cent. while Fig. 5 gives the power obtained from such 
a windmill with a chord of 2 cm. as approximately 0.54 kilowatts. Now Fig. 7 
shows that the maximum efficiency of an airscrew and the power obtained from it 
increase as the angular velocity increases. This suggests that, although the maximum 
stresses will also be increased, higher efficiencies can be obtained by increasing the 
angular velocity. In any particular case, it would clearly be desirable to work out 
all these results for a range of values of 2c. The method used for integrating the 
equations is a very general one so that the theory can be applied to a wide range 
of airscrews and, in particular, the straight leading edge and constant cross section 
chosen for simplicity in the example worked out in this paper are by no means 
essential. Accounts of the methods used in propeller design which presents problems 
similar to those in the case considered in this paper may be found in such works as 
those by Theodorsen“”’ and Bairstow"). 
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In cases where greater accuracy is required than was obtained in the particular 


example considered, it is always possible to calculate further terms of the series 
expansion for the potential in powers of x. Terms in higher powers of x can then be 
retained in calculating the torque and the drag and the efficiency and the power 
obtained to any required degree of accuracy. 
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Pohlhausen’s Method for Three-Dimensional 
Laminar Boundary Layers 


J. C. COOKE, M.A. 
(Professor of Applied Mathematics, University of Malaya) 


(Summary: Wild has extended to three dimensions the von Ka4rman-Pohlhausen 
method for the laminar boundary layer flow over a fixed obstacle and used the 
method for an infinite yawed elliptic cylinder in a stream. In this paper the method 
is tested in two ways (which may be called full-Pohlhausen and semi-Pohlhausen) for 
the case of an infinite yawed cylinder when the velocity outside the boundary layer 
over the surface normal to the generators is of the form U=cx". The exact 
solution is known in this case. 


A table of the skin friction, displacement thickness and momentum thickness 
is given for various values of 8 [=2m/(m+1)], and the agreement is found to be 
fairly good for 8>0 (accelerated flow) but not so good for 8<0 (retarded flow).) 


1. Introduction 


Wild” has given a method for calculating the three-dimensional laminar 
boundary layer over an infinite yawed elliptic cylinder by an extension of the method 
of Pohlhausen and von Karman. While this method is satisfactory under a favour- 
able pressure gradient, it is not so good when the flow is retarded. Nevertheless it 
gives useful results and so the extension of the method to yawed infinite cylinders is 
worth consideration. This paper considers yawed infinite cylinders when the 
potential flow normal to the generators is of the form U=cx”. The exact solutions 
for this case are known’? and so it is possible to check the method. The results 
are fair for positive values of m, less good for negative values. Even in the simple 
case given by Wild the numerical work is fairly heavy, although straightforward. 
It would seem that it will be very heavy indeed if the method is extended 
to more general cases. 


Notation 
x,y,z rectangular Cartesian co-ordinates, the x-axis being tangential to the 
surface and normal to the generators, the y-axis being parallel with 
the generators 


Paper received January 1951. 
[The Aeronautical Quarterly, Vol. III, May 1951] 
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u,v,w velocity components in the boundary layer in the directions of 
x, y and z respectively 


vv first and second derivatives of v with respect to z 


U,V main stream velocity components in the directions of x and y 
respectively 


U’,U” _ the first and second derivations of U with respect to x 
c,m constants in the equation for the velocity U=cx” 
B 2m/(m+1) 
p density 
Pp __spressure 
coefficient of viscosity 
vy kinematic viscosity = u/p 
8, boundary layer thickness of the velocity component v 
6, boundary layer thickness of the velocity component u 
A 
f(A), g(A) functions defined by equation (7) 


8.* displacement thickness of boundary layer for velocity component u 


(= | (-u/U)dz) 
0 
5,* displacement thickness of boundary layer for velocity component v 
by 


(= | 
0 
2 7y components of skin friction in directions of x and y respectively 
 4,/8: 
4(m+1) 
T, 


» momentum thickness of boundary layer for velocity component v 


0 
4% k/[U/(x)] 
D* 8,* /[+mU /(vx)] 
= 8,* ki [U/(»x)] 
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(U/(»x)] 

Y/D 

integral defined immediately after equation (10) 

a function of K and 2 defined by the relation 


d 
> 
ax (FU 6,) = —2¥/(K 6,) 


2. The Equations of Motion 


Taking the x-axis tangential to the surface and normal to the generators, the 
y-axis parallel to them and the z-axis perpendicular to those of x and y, the Prandtl 
equations of boundary layer flow become 


dv/dy+ dw/ dz 


The main stream velocity components are taken as U and V, where V 
is constant. 


In the case of infinite yawed cylinders none of the functions changes as y 
changes, and we have 


(1), (3) and (4) are the ordinary boundary layer equations for two-dimensional 
flow in the xz-plane, of which there are many methods of solution including 
Pohlhausen’s. (2) is a new equation for v, the cross flow in the boundary layer. 


53 


of D 
Y 
a 
Ou Ou Ou lop 07u 
> 
Ox cy pox 02? 
ov ov ov ldp 
Ox Oy p Oy 
u 
O=dp/dz, 
0. 
y 
Ox pox 
ov ov 
Ox 0z 
4 


J. COOKE 


By integration of (2) with respect to z, or directly from momentum considerations 
(see also Tetervin’) we obtain as a momentum equation 


| 


~) 
0 


The von Karman-Pohlhausen assumptions are (a) that this equation holds over 
a small boundary layer thickness 5,, where here 3, is a function of x only, (b) the 
velocity v may be expanded in powers of z/5, up to the fourth power. When the 
usual boundary conditions y=v”=0 when z=0; v’=v”=0, v=V when z=4, are 
applied (the quartic will not of course fulfil all the boundary conditions) we obtain 


v/V =2(z/8,)—2 + (z/8,)* 
as given by Wild. 


The two-dimensional Pohlhausen method is too well known to be described in 
detail here. The best account is given by Howarth’ who repeated it in a shortened 
form in Ref. 5. Suffice it to say that u/U is expanded as a quartic in z/é,; the 
boundary conditions lead to the equation 


+5) 
where A=U' 3,?/v. 
Substituting in the Karman momentum equation leads to 


dashes denoting differentiation with respect to x. f(A) and g(A) are known rational 
functions of A. 


Here the velocity external to the layer is considered as U=cx", which makes 
(7) become 


dd /dx=(1/x) { mf (A)+(m—1)gQ)}. 


This equation can be integrated directly and it appears that it has no solution 
except A=constant; this is determined by solving the equation 


mf (A)+ (m—1)g(A)=0. (8) 


That A should be a constant is expected by similarity considerations; the 
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solution is exactly known (Falkner and Skan‘*’”) and Hartree‘*)) and the velocity 
profiles are known to be “ similar.” 


Equation (8) is a cubic and the correct root is determined by physical 
considerations in more general cases. Dryden" gives a curve of A against 8, 
where B=2m/(m + 1). 


It is worth while showing the accuracy of this solution before proceeding to the 
three-dimensional case, and the results for the skin friction +, and displacement 
thickness 4,* are given in Table I for selected values of 8. In this table 


D.*=8,* /[+mU/(x)]. 


It can be seen that the accuracy is reasonable except for the skin friction in the 
case of negative 8. It is well known that this method gives inaccurate results in a 
retarded region. 


3. The Equation for v 


by 


We have ax 


where v is given by the equation 


for 0 <z< 4, together with y=V for z > 4,. 


Also u/U is given by the quartic (6) in z/é,, if 0 < z < 4, together 
with u=U if z > 4,. 


TABLE I 


D.* 


(approx.) (exact) (approx.) (exact) 


0.326 0.000 0.936 1.052 
0.886 1.010 0.469 0.446 
1.527 1.536 0.442 0.440 
1.254 1.286 0.588 0.592 
1.196 1.233 0.641 0.648 
1.164 1.202 0.679 0.688 
1.158 1.193 0.693 0.703 


e 
n 
d 
S B A 
—0.1988 — 6.869 
~0.1 ~2.170 
0.2 2.509 
0.6 5.268 
. 1.0 7.052 
1.6 9.349 
e § 2.0 11.984 
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The integral will be different according as 5, is greater or less than 3, owing 
to the differences in the forms of u/U and v/V inside and outside the boundary 
layer. u/U and v/V do not fulfil all the boundary conditions. 


Writing 6,=K 45, we obtain after some reduction 
d 2v 
dx K 5,’ 
where F is a function of K and A. 


This function has the value 


F, =0.3 — 0.00833 A + (1 / K) (— 0.13333 + 0.00556 A) + (1/ K*) (0.02143 — 0.00119 A) + 
+(1/K*)(— 0.00556 + 0.00033 0.3K 


if K>1, and the value 


F,= — K? (0.13333 0.01111 A) + K* (0.01190 A)+ K* (0.02143 — 0.00536 A) + 
+ K* (— 0.00556 + 0.00093 A) 
if K<1l. 
These results agree with those of Wild, although written differently. 


Using the equations 6,?/v=constant, and U=cx”™ we obtain 
dK |dx= —(2/K+AF/8)](xdF/dK). 


At the forward stagnation line x=0 and if dK /dx=00 is excluded as physically 
impossible, the numerator must also vanish for x=0. The equation has in fact a 
singular point at x=0, and in all the instances tried it is found that the only solution 
of the equation is K=constant, the constant being determined from the equation 


28+AKF=0. 


In the case K>1 this is a quintic for K and in the case K<1 it is a sextic. In 
each of the cases tried there is only one physically acceptable root in the 
desired range. 


4. Results 


The results have been worked out for the values —0.1988, —0.1, 0.2, 1.0, 
1.6, and 2.0 for 8, and are given in Table II. To obtain results for the skin friction, 
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and writing 


For the displacement thickness 
by 


0 


For the momentum thickness 


y 


0 


B K 
(approx.) 
—0.1988 —6.869 0.8574 0.3968 
—0.1988 (Alternative method) 0.360 


—0.1 —2.170 0.9412 0.4561 

—0.1 (Alternative method) 0.460 
0.2 2.509 1.089 0.5188 
0.6 5.268 1.2154 0.5555 
1.0 7.052 1.3043 0.5773 
1.6 9.349 1.3819 0.5985 


2.0 11.984 1.3400 0.6099 


THREE-DIMENSIONAL 


direct comparison can be made with 


writing D,* =5,* kV (U/(vx)] 


0 


TABLE II 


T, 
(exact) 
0.3258 
0.3258 
0.4332 
0.4332 
0.5069 
0.5467 
0.5704 
0.5938 
0.6053 


BOUNDARY 


a direct comparison is possible with G’(0) in Ref. 2. 


a¢= | 


6, = dz: 0.11746 8,, 


D,* 
(approx.) 
1.512 
1.667 
1.316 
1.304 
1.1565 
1.0801 
1.0393 
1.0025 
0.9828 


LAYERS 


0.3K /(A/B), 


(1 —G)dY in Ref. 2. 


and writing 0,=6,k/[(U /(vx)]=0.11746 K /(A/8), 


a direct comparison can be made with | G(1—G)dY in Reef. 2. 


(exact) 


1.6696 
1.6696 
1.3043 
1.3043 
1.1376 
1.0646 
1.0265 
0.9917 
0.9759 


[k?=4(m+1)] 


0, 
(approx.) 
0.5920 

0.653 

0.5160 
0.511 

0.4497 
0.4173 
0.4044 
0.3923 
0.3872 


0, 
(exact) 
0.5937 
0.5937 
0.5151 
0.5151 
0.4528 
0.4229 
0.4069 
0.3925 
0.3852 
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It is seen that the agreement is quite good for positive 8, that is for accelerated 
flow in the x direction. Even in the case of retarded flow the agreement is better 
than in the two-dimensional case. As it depends on the two-dimensional solution, 
which is so inaccurate here, it cannot be expected to be good. 


5. Alternative Procedure 


In view of the foregoing it might seem better to use the known two-dimensional 
solution, and to use Pohlhausen’s method for the cross flow only. 


In the equation 


a v-V)de= 
u(v—V)dz=—y 


0 


the Pohlhausen quartic (9) is substituted for v/V and the Hartree value UF’ (Y) for 
u, where Y =k /[U/(vx)] z, k?=4(m+ 1) (the notation is that of Goldstein’). 


On putting D=k/[U/(vx)]5,, £=Y/D, and assuming in virtue of the known 
“similarity ” that D is constant we obtain after some reduction 


1 
where I= ) (-26+26°- DF’ (DOdé. 


0 


To solve (10) D is given a set of suitable values, 7 is evaluated numerically 
using Hartree’s table, and by interpolation the value which satisfies equation (10) 
is found. 


In this paper the method has only been used approximately for the values 
—0.1 and —0.1988 for 8. It is found that 


The values are given in Table II. It will be seen that while some of the values 
are improved others are not so good as before. As can be seen from equations (11), 
the three quantities tabulated depend on each other and it is impossible to make 
them all agree at the same time. It shows in fact that the fundamental assumption 
of a quartic is not quite good enough in this region. 


If the velocity profiles are plotted it is found that the curves for the two methods 
(exact and approximate) run very close together for all the values of 6 tried here 
except 6= —0.1988. In Fig. 1 the curves show the exact solution, the first method 
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v/V 4 
0 4b 
VELOCITY PROFILE FOR = 01988 
: EXACT 
FULL- POHLHAUSEN 
—-— SEMI-POHLHAUSEN 
k /[u/(x)]z 
Fig. 1. 


given here (full-Pohlhausen) and the alternative procedure (semi-Pohlhausen) all 
for 8=—0.1988. These curves show clearly that the assumed quartic velocity 
profile is not sufficiently good in the region of retarded flow in the main stream. 


6. Conclusions 


It seems doubtful if the accuracy of the method would be worth the labour 
involved in more general cases. It is possible theoretically, using the momentum 
equations of Tetervin’, but the labour would seem to be prohibitive. Prandtl”° 
has suggested a method, but these methods would require testing by knowledge 
of more “exact” solutions than are at present known. Apart from flows involving 
axial symmetry, and from yawed infinite cylinders, there appears to be only one 
exact three-dimensional solution known. (Carrier 


Taylor’*’ and Binnie and Harris"*) have used Pohlhausen’s method for the 
solution of swirl problems but it may be remarked that they assume that lengths 
equivalent to what are called 5, and 5, in the present paper are to be taken as equal. 
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It is certainly true that in the exact solution given in Ref. 2 the boundary layer 
thicknesses in the two directions (reckoned on any consistent basis, e.g. the usual 
99 per cent. of the main stream velocities) are considerably different. Their ratio 
in fact agrees fairly well with K of this paper. It is hoped to discuss this matter in 
a later paper. 
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Supersonic Flow Past Bodies of Revolution 
with Thin Wings of Small Aspect Ratio 


P. M. STOCKER 


(Department of Mathematics, The University, Manchester) 


(Summary: The method developed by G. N. Ward" for the treatment of slender 
pointed bodies in a uniform supersonic stream is applied to three special cases. 


({) Supersonic flow past a body of revolution with thin wings of symmetrical 
section and of small aspect ratio at zero incidence. 


(ii) Supersonic flow past a body of revolution with plane wings of small aspect 
ratio set at incidence to the body, the whole being at incidence to the stream. 


(iii) Supersonic flow past a body of revolution with a plane fin of small aspect 
ratio set at incidence, the whole being at incidence to the stream. 


The pressure distribution on the wing has been calculated for a special case of (i) 
and is given in the Appendix.) 


1. Introduction 


In Ref. 1 the well-known linearised equation for the supersonic isentropic flow 
of a perfect inviscid compressible fluid past slender pointed bodies of thickness 
ratio t (t being small compared with unity) is solved by operational methods. It is 
shown that if the perturbation potential is U¢, the function ¢ may be approximated 
to within a factor 1 +O(f log t) by the function ¢, obtained as follows. Taking any 
section of the body by a plane perpendicular to the undisturbed stream, a solution 
¢, of the two-dimensional Laplace equation V? #,=0 is determined, which satisfies 
certain boundary conditions on the boundaries of this body section. The boundary 
conditions which must be applied are such that the solution is found to be identical 
with the solution of a problem in two-dimensional incompressible flow. If the 
plane of section is moved downstream with velocity U, the shape and position of 
the body section relative to the plane will change; this may be regarded as the 
motion of a cylinder of varying cross-section in an incompressible fluid at rest at 
infinity. The potential 9, at any point is equal to the potential at that point, in the 
plane of section through the point, due to the incompressible motion just described. 


Paper received June 1950. 
[The Aeronautical Quarterly, Vol. III, May 1951] 
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A system of right-handed orthogonal axes Ox, Oy and Os is chosen with the 
origin at the nose of the body, the axis of s being parallel to the undisturbed stream 
and extending downstream; the scale of s is chosen so as to make the body of 
unit length. Let S(s) be the area of the section of the body by a plane s=constant, 
let z=x+iy and let w,=9¢,+i,, where ¥, is the harmonic conjugate of 9,. For 
sufficiently large |z|, w, may be expanded in the form 


=) =a, logz+b, +3 ay 
=1 


where a,, b, and ay are functions of s. Let x,(s) and y,(s) be the co-ordinates of 
the centre of area of the body section in the plane s=constant. Let p, be the density 
of the fluid in the undisturbed stream and z,=x,+iy,. It is shown in Ref. 1 that if 
X and Y are the components of total force on the body in the directions of x and y 
increasing respectively, then 


F/(4p, U*)=(X +iY)/G p, (a,),-1 +28’ (1) z, (1) +28 (1) z,’ (1), (1) 


where primes denote differentiation with respect to s. 


The drag force D in the direction of s increasing is given by 


log —— (s) S” (c) do ds— og S” (7) do — 
0 


2n 


where c,x35(1) is the base drag. The pressure p is given, to a consistent order 
of approximation, by 


where z=re”. 


Strictly, this theory can only be applied to bodies of maximum thickness 1, 
where ¢ is small compared with unity, for which the angle which any tangent plane 
to the body boundary makes with the undisturbed stream direction is small and 
O(t) and for which the rate of change of this angle is O(t). Moreover, if d is the 
maximum diameter of a section of the body by a plane perpendicular to the stream, 
then the maximum curvature of this section must be O (d~') for all points where the 
section is convex outwards. This last condition is not met by the winged bodies 
considered in this paper, but it is shown in Ref. 1, p. 90, that the solution obtained 
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is an approximation to the true potential except in the neighbourhood of points 
where the curvature is infinite (i.e. the wing edges); furthermore, as the pressure 
singularities are integrable the expressions for the body forces should be true 
approximations for small incidences. It is difficult to say for what range of aspect 
ratio the theory is valid and it would seem that this can only be determined by 
means of an extended theory or experimental evidence. 


Notation 
U __ velocity of the undisturbed stream 
t _ thickness ratio of the body 
perturbation potential 
U¢, approximate perturbation potential 
v¥,v¥, complex conjugates of and ¢, 
Wy 
x,y,S rectangular Cartesian co-ordinates 
z x+iy 
a,,b,,@x functions of s only (N from 1 to 0) in the expansion of w, (see 
Introduction) 
S(s) cross-section area of the body and wings by a plane s=constant 
Xg.y,_ co-ordinates of the centre of area of this section 
Ze Xgtiy, 
I‘ contour of this section in the plane s=1 
v,7 Outward normal and tangent to I respectively 
p, density of the fluid in the undisturbed stream 


X,Y components of total force in the directions of x and y increasing 
respectively 


F X+iY 

D drag force in the direction of s increasing 
coefficient of base drag 
o,X variables of integration 

p pressure 

Pp, pressure in the undisturbed stream 


r,6 polar co-ordinates in x, y-plane 


d maximum diameter of body section by a plane perpendicular to 
the stream 
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thickness of the wing at (y, s) 
cross-section area of body alone 
body radius 

wing semi-span (Sections 2 and 3); height of fin (Section 4) 
outward normal to body and wings 

— (Sx’ (s) log r)/(2™) 

strength of line source along x=0 


: > (see Section 2) 
strength of image source j 


real and imaginary parts 

Mach number of undisturbed flow 

Z+as 

real and imaginary parts of z, 

incidence of body 

incidence of wings (Section 3), of fin (Section 4) 
B—a, incidence of wing to body 

complex potential corresponding to ¢, 

W,—-W, 

polar co-ordinates in x,, y,-plane 

value of corresponding to w, 

real and imaginary parts of ¢ 

2c 

a+c*/a 

09,/0€, 0¢./0n 

functions of R’, a’, z and 6 (see equations (18) and (19) ) 
lift of system 

lift when 

extra lift due to incidence y of wing (or fin) to body 
lift of body alone 

lift of wing alone 

lift of fin alone 


drag on system in position of zero lift 
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2’ incidence of system from zero lift position 
D, drag on body alone 
2, incidence of body in zero lift position 
D, drag when B=0 
6b maximum value of a(s) 
R, _ value of c(s) where maximum value of a(s) occurs 
K,,K, functions of b and R,, defined between equations (23) and (24) 
B,,C,,d,,d,,d, functions of z, 8, R’ and a@ defined by equations (28) and (29) 
K constant of integration in equation (30) 
S.(s) area of wing section by plane s=constant 
Fig. 2 and Appendix 
k (1/=)(da/ds) tan 
a*,y* ale, y/e 
bos. yw contributions to b, from body and wing respectively 


Px _— pressure due to body alone 


2. A Body of Revolution with Wings of Symmetrical Section and 
of Small Aspect Ratio at Zero Incidence 


The wings are assumed to be symmetrical about their chord plane x=0 and 
to be sufficiently thin to enable the boundary conditions on the wing surfaces to 
be applied on x=0 instead, i.e. their thickness must be small compared with the 
body radius. Let the radius of the body be c(s), the semi-span of the wings a(s), 
the thickness of the wing at the point (y, s) be 2f(y, s) and the area of the body 
section Sy (s). 


The appropriate two-dimensional problem is as follows. To find a solution 
¢, of V?@,=0 in the plane s=constant such that 


_ s)]=f' say, on c<|y|<a,x=0 (5) 
on = y, = y,S y, é 


where n is the outward normal to the boundary of the section of the body and 
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wings by the plane s=constant. The dependence of a, c and f(y) upon s does not 
affect the solution of the two-dimensional problem and indication of the functional 
dependence will be omitted. The boundary conditions may be satisfied by placing 
a distribution of sources and sinks along the line x=0. To satisfy condition (4) 
a source of strength S’; (s)/(2) is placed at the origin. If 9, —(S’s (s) logr)/(2z), 
the function ¢, satisfies the boundary conditions; 


(y) on c<|y|<a, x=0, (8) 
0o,/dr—>0 as r—>oo. ‘ ‘ (9) 


¢, is determined by placing a line source of strength F (y) along the line x=0 from 
y=-—ato y=-—c and from y= c to y=a. To satisfy condition (7) an image line 
source must be added from y= —c to y= —c*/a and from y=c’/a to y=c together 
with a sink of appropriate strength at the origin. The strength, ®, of the image 
source is given by 


F (y)dy= —®(c*/y)d(c*/y), 


ie. F()=(c?/y*) ®(c*/y), 
and the strength of the sink by 


a 


| F (y)dy. 


c 


F (y) is now determined in terms of f(y) by the use of equation (8). Let w, be 
the complex potential corresponding to ¢,; then 


c 


| Flog (22 +A*)da+ dd —2 log z 


c*/a c (10) 
dw, F(A) c’ F (A) 2/ 


Equation (8) gives 


R =f (y) on c<|y|<a, x=0+. 
Z Iz=iv+0 
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From (11), assuming that c<|y|<a for definiteness 


= ict y | + = | F (A) dA 2iy | + Fy), 
J y F ( )d + (A?— y?) F (y) 


where the Cauchy Principle Value is to be taken in the third integral. 


Thus the complex potential w, is given by 


2 jogz | f (s)logz+by. (13) 


The term b, is a constant in the two-dimensional problem, but is in general a function 
of s. It is shown in Ref. 1 that 5, is given by 


Ss 


b,= { tog (s—o)S” (v) do } » « 


0 


where M is the Mach number of the undisturbed flow. 


Knowing w, and b, the pressure distribution on the body and wings can be 
calculated, and by a combination of this solution with that obtained in Section 3 
for plane wings at small incidence the pressure distribution is obtained for a 
body of revolution with wings of symmetrical section at small incidence. 


The pressure distribution on the wing has been computed for a body of 
revolution with wings of parabolic plan form and the results are given in 
the Appendix. 


3. A Body of Revolution with Plane Wings of Small Aspect Ratio 
at Incidence, the Whole Being at Incidence to the Stream 


Let = be the incidence of the body to the stream and £ that of the wings; 
it is assumed that |8—«| is sufficiently small to permit the boundary conditions on 
the wing surfaces to be applied on the plane x= — as instead. It will be necessary 
to assume, as in Ref. 1, that the body is of one of the following forms: 
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(i) a body of revolution of general variation of cross-section with wings of 
which the edges are everywhere leading edges, upstream from the base; 


(ii) a body of revolution of uniform cross-section downstream from the section 
of maximum wing span. 


These restrictions are necessary because downstream of the section of maximum 
span a vortex sheet is formed and the determination of the solution in the general 
case when the body radius is not constant there is very tedious. 


Let z,=x,+iy,=z+as and let 6,=argz,. The boundary conditions on w, 
in a plane s=constant which meets the leading edge of the wing are 


09, /dn= —acos 6,+S’ (s)/(2xc) on |z,|=c, 
on x,=0+, c<ly,|<a, (15) 
l\dw,/dz,|=O(1/|z,|) at infinity . 
The case 2=8 is solved in Ref. 1, the solution being 
w, — b,=[1/(2)] (s)log z, +8 { [(z, +¢7/z,? } . 
If w,—w,=w.,, the boundary conditions on w, are 
09,/dn= —(a—f)cos 4, on |z,|=c, 
on x,=0, c<|y,|<a, : . (16) 
ldw,/dz,|=O(1/|z,|*) at infinity. 
The conformal transformation 
(=£+in=z,-c?/z, 


is now applied. 
The body section is mapped onto the imaginary axis between 


(= —(a+c?/a)i and (=(a+c*/a)i. 
|\d&/dz,|=2 |cos on |z,|=c. 
The boundary conditions on w, in the ¢-plane are therefore 
on £=0+, |n|<2c, 
0¢,/0n=0 on £=0, 2c<|n|<a+c’/a, (17) 
at infinity. 
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Let R’=2c, a =a+c*/a, u,=09,/0£, v,=09,/0n and let g,=u,—iv, be the complex 
velocity. The boundary conditions (17) show that du,/0n=0 on €=0, |n|<a@’ and 
therefore dq,/d¢ must be real on this segment. gq, would be expected to have a 
branch point at the tips of the wing and dq,/d¢ to have poles at the points where 


the wings meet the body. Hence, remembering that |dq,/d¢|=O(|¢-?}), 


dq. A’ 


(18) 


where A’ and B’ are real constants; it can be seen from symmetry that there is no 
term in ¢ in the numerator. 


Integrating, we obtain 


a 


q2>= (+ ay +R? 
A is real and B is imaginary. A, B and C are now determined from the boundary 
conditions. The jump in the imaginary part of the term cosh~'[R’((?+a’)}/ 
| a’ (? +R”) } ] on the upper surface of the wing at (=iR’ is 477. Thus, since the 
jump in & [q.] is 4(a—8), 
B=i(8—2)/z, 
and, since #[g,]=4(8— 2) at (=0, 
A(C]=4 (8-2). 
The condition g,=0 at infinity yields 
A+Bcosh-' (R’/a’)+C=0, 
that is A=(2—) {4—(1/) cos~' (R’/a’)} and .A[C]=0. 
It follows that q, is given by 
q.=(B-2) { (Leos (R’/a’)—4 — cosh } 
(20) 


To find », further integration is necessary, but the lift can be obtained without 
integration since only a,, the coefficient of 1/z in the expansion of w, need be 
known. The coefficient of (-? in the expansion of g, in descending powers of ( is 


A+4iR’ B(a’* — 
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which is equal to 
{ (8—2)/(2=) }{ sin-' (R’/a’)— R’ (@?- } . 


This coefficient is equal to the coefficient of —1/z in the expansion of w,. The 
coefficient of 1/z in the expansion of w, is 


(a +c*)/a’, 


and hence 


(21) 


As mentioned earlier, this expression holds only so long as the wing edges are 
leading edges; however, it is shown in Ref. 1 that the solution is of the same form 
downstream of the section of maximum span, but a, the semi-span, must be replaced 
by b, the maximum semi-span and the body radius c by its value, R,, at the section 
of maximum span. The same is true in this case; the strength of the vortex sheet 
is independent of s and hence of a. Thus a solution is obtained which is identical 
with that for a wing which has constant span downstream of the section of maximum 
span. It should be noted that this is only true so long as one of the conditions 
(i) and (ii) given at the beginning of this Section is met. 


The lift on the system is given by 


L 
+ 2S’ (1) z, (1)+ 2S (1) z,’ (1) 
(22) 
This may be written 
(b7+R,?¥ . _,f 2bR, 
4p,U* 4p, U? b? 
b*+R,* 


where A is the lift when both body and wing are at incidence «. Thus the extra 
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lift, L,, due to the incidence y of the wing to the body is given by 


The effective lift of the wings is always increased by the presence of the body, 
so long as the incidence of the wing is greater than that of the body. If L, Ls, 
and Lw denote the lifts of the system, body, and wing respectively, we have"? 


Ly/(4 p, U*)=22 R,? 2, p, B. 


From (22), it can be seen that 


L-L, E —2R,? +R,‘ 4 -1( 2bR, 
B ] +2(2—f)| — sin 


(b*= 


Hence, using the fact that z< and the inequality sin~' x < 4 =x, we have 
L-Ly > Lw(1+R,/b+R,?/b’). 
The draz may be calculated from the formula (Ref. 1, p. 90) 
D=D,+342'L, 


where D, is the drag in the position of zero lift and 2’ is the incidence of the system 
from that position. Writing equation (22) as 


L/(Qp, U*)=K, a+K,y, 


where K,=22 (b*-—b? R,?+R,*)/b? 


. ( 2bR, ‘ 
and K, sin += (b?-R, )]. 


we have, if 2, is the incidence of the body in the position of zero lift 


K, a,+K, y=0. 
Hence L/G Pi U*)=K, 2)=K, a’, 
and thus D/(4p, U?)=D,/G p, U?)+4 K, . . (24) 
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To determine D, it is most convenient to determine first D,, the drag when 
B=0. From equation (2) it may be seen that D, is the sum of two terms, the drag 


of the body alone at zero incidence and that part of the integral — | oo (0¢,/0v)dr 


which depends on a. In general the wing (or vortex sheet) must be included in the 
contour I"; however, in this case the integral over the wing is zero, since on it 
0¢,/0v=+B8=0. The contour I may be taken to be the circle z,=R,e”; on 
this circle 


d>=R, d6,= —iR, dz,/2z, 
and 0¢,/dv= —2cos 0, +S’ (1)/(2= Rj) 


= — 2 (2,7 +R,*)/(4z, R)+S R,). 


tus 9, 4 Jim, de] 
r 


2 


The term dependent on z is found to be 
—%a(a,)e-0+4 27S’ (1). 
Hence, D,/ p, (@,)g-0—4 2 S (1) + p, U’), 


where Dy, is the drag on the body alone at zero incidence. Using equation (24), 


4p,U? 4p,U? "2 b* b?+R,? b K, 
Pr 1 


Therefore for any given body and given wing, if the lift is prescribed, the drag is 
a minimum when the body and wing are at the same incidence to the stream. 
To see this it may be noted that a’=L/K,, and hence for given lift «’ is constant 
and from equation (24) the drag has a minimum when D, has a minimum, that is, 
when y=0. 


4. A Body of Revolution with a Plane Fin (i.e. a Wing on One 
Side Only) of Small Aspect Ratio Set at Incidence, the Whole 
Being at Incidence to the Stream 


The method of solution is identical with that of Section 3. Let the incidence 
of the fin to the undisturbed stream be 8 and that of the body 2; let the body radius 


72 


= 


SUPERSONIC FLOW 


be c(s) and the height of the fin a(s). Let the fin lie in the plane x= — as and let 
z,=xX,+iy,=z+2s. The boundary conditions in the appropriate two-dimensional 
problem are 


0¢,/On= FB on x,=0+, c<y<a, 


29, —acos (s)/(2xc) on z,=ce®, | 
(26) 
J 


/dz,|=O(|z,|~") at infinity . 


The case B=a is first solved. Let the solution to this problem be 
w,-22,+S (s)logz,/(2™); w, satisfies the boundary conditions 


dw,/dz,= at infinity , 
on z,=ce®, and on x,=0,c<y, Sa. 


If the transformation (=z,—c*/z, is now applied the problem is reduced to that 
of the flow of a uniform stream of velocity « past a flat plate at right angles to the 
stream; the length of the plate is a(1+c/a) and its mid-point is at (=4ia(1—c/ay’. 


Thus w,=a[ { } 

As in Section 3 this holds only so long as the edge of the fin is a leading edge. 
To find the solution downstream of the section at which the maximum height of 
the fin occurs the same assumptions must be made regarding the body section as 


in Section 3; the solution is then found by replacing a by its maximum value b, and 
c by its value R,, at the section at which a attains its maximum. 


The coefficient of 1/z in the expansion of w, when s=1 is 


b? 
a[ (14% —R, 


The coefficient of 1/z in the expansion of the full potential in descending powers 


of z is 
b* R,\* 2 


and hence the lift is given by 


L/ (4p, U2)=4e2 { $b? } . 
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The lift of the body alone is given by 
Ls/(4 p, 
and that of the fin alone by 
Ly/(4 p, U*)=4 22 
Thus 
Returning to the case when 82, it is sufficient to solve the problem 
on |z,|=c, 
on x,=0,c<y<a, 
|\dw,/dz,|=O(\z,|-?) at infinity, 


for by adding the solution of this problem to the one obtained for the case B= 
the general solution is obtained. 


Again applying the conformal transformation 
the bounc.ury conditions become 
0o,/On= +4(B—2) on £=04,-—R’<n<R’, 
éo,/dn=0 on £=0, R’<n<ad, 
ldw,/d{| =O(|{|-*) at infinity, 


Let gq. be the complex velocity dw,/dt. q, is expected to have branch points at 
(= —iR’ and (=iad and dg,/d{ to have a pole at (=iR’. Hence we write 


dq, 


28) 
de ~ 
On €=0, —R’<n<d’, 0q,/0n is imaginary and hence B, and C, are real. 
d d d 


d,, d, and d, are determined from the boundary conditions and from the equation 


d, +d,+d,=0. 
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Integrating dq,/d¢ it is found that 


d, lo 3a’ R’ — R” [2R’ (a — R’ _ 
2d, iv 2d, iv 


+K, 


~ 
(30) 


where K is a constant of integration and YW =[((+iR’)((—ia’)]}!. The jump in the 
real part of g, at (=iR’ is 4 (z—{) and this gives 


d, =[(2—)/(2=)] [2R’ — R’)}*i. 


At(=0, A[qg.]=4 (8-2) and hence 


[25 itog 3a’ R’ — R® — 2R’ + FAK] =4(8-2), 
2m —iR’ 

i.e. AR{K]=} (8-2) 

The boundary condition at infinity yields 


{log + GER) 


+K=0. 
Taking real and imaginary parts 

and /[K]= log (a’+R’). 


The coefficient of 1/z in the expansion of w, is 


{ arg{ RY} | +=/2}. 


Thus the extra lift due to the extra incidence, 6 — 2, of the fin to the body is 


[are({°*: —4 E (6-R, |’) +2/2]} 
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APPENDIX 
The pressure distribution on the wings of a system consisting of a body of 


revolution with wings of symmetrical section, at zero incidence, has been calculated 
in two special cases. 


The dimensions of the body are taken as follows. The cross-sectional area of 
the body is defined by 


Sz (s)=(0.05)? = + (0.08) = (s—0.5)[1 + 1.5 (s—0.5)], if O< s <0.5 
=(0.05 =, if 05S<s5 <1. 


The body is chosen so that S$, (s), S’3(s) and S$”, (s) are continuous at s=0.5. The 
wing span, a(s), is defined by 


a-—c=0, ifs=<0.5 
=1.6(s—0.5)(l-s), if <1. 
The section of a wing by a plane s=constant is taken to be a triangle, with angle 


2u at the tip, and « is taken as tan-' (0.1). The plan form is shown in Fig. 1 and 
a section in Fig. 2. 


0-1 

0:1 

| 
| | | 
0-5 +" 0-5 

Fig. 1. 


é 
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Fig. 2 


From equation (13) we find that 


Ob, _ #2_ ye2)_ 9 * a* 
[a log (a** — y**) — 2 log (y** — 1)— y* log| 


+y*log(% 
2(a*-1 * #2 #2 
—2 (a* —1)—2(a* — 1) log y+ log + a* log (a** y** 1)— 


d 
2a log a* +2 (a*— 1) loge | +k Flog (a*? y**)+ 
+log (a*? y*? — 1) — 2 log y* —2loga* +2 ogc], 


where k=(1/%)(da/ds)tanu, a*=a/c and y*=y/c. 


and 09,/0x=7k. 6, is given by formula (14) and denoting by b,w and bd,» the 
contributions to b, from the wings and body respectively, we find that 


5 


b’.3= (s) log [4 (M? — log (s—o)S”’s (c) do —S”, (0) log s 
0 
and 


(s) log [4 (M? — 1)']— | log (s—o) do — (0.5) log s } 


0.5 


In the computation the Mach number of the undisturbed stream is taken to be 
/5 and hence (M?—1)'=2. 
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The results are shown in Figs. 3 and 4. Lines of constant pressure on the wing 
are shown in Fig. 3, the values of (p—p,)/(4 p, U*) being indicated. In Fig. 4 lines 
upon which (p— ps)/(4 p, U*) is constant are shown, ps, being the pressure due to 
the body alone. The solution gives an infinite pressure at the wing edge, and the 
isobars have not been computed in the neighbourhood of the edge. 
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Normal Elliptic Functions: A Normalized Form of Weierstrass’s 
Elliptic Functions by A. R. Low, M.A. University of Toronto 
Press. Toronto, 1950. Price $1.25 (Can.). 


The general quartic, with four complex roots, may be transformed by 
homographic transformation to a cubic with three complex roots and one root 
at infinity. The further arbitrary condition that the sum of the roots is zero gives 
Weierstrass’s standard form. When the roots of the cubic are real the form is that 
which arises naturally in a class of physical problems referred to a system of 
confocal real conicoids, ellipses, one-sheet hyperboloids and two-sheet hyperboloids 


The quartic was normalised by Legendre to the form 
S=(1 —s*) (1 —k? s*) <4. 
with only one real arbitrary parameter k. 
In analogous manner the cubic is normalised by the author to the form 
P=p(p—m)(p-—1) 0<m<l. 


Taking m=k?, the substitution p=s~? transforms Weierstrass’ elliptic integral 
with normalised cubic to Legendre’s first elliptic integral 


u=3| P-tdp= | s ‘ds. 
0 


The real p-axis is divided into four segments by the roots of the normal cubic, 
0, m, 1, and the concealed root at infinity, and four elliptic integrals are defined, 
one in each range, denoted by suffixes, p,, p., Ps, P,, being regarded as independent 
variables. Inverting the functional relations between u and p, four elliptic functions 
are defined and may be written p, (u,m), p, (u,m), p, (u,m), p,(u,m), to denote 
the functional nature of p explicitly, u being now regarded as the independent 
variable. 


The equation of the confocal system becomes 
x? p-'+y? (p— 
Putting p=1, a limiting form is obtained as an elliptic lamina 


x*+y?(l—m)-'=1 
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which may be taken as the conicoid of reference, with unit semi-diameter and only 
one arbitrary parameter, m=eccentricity squared. 


This is analogous to referring the elementary transcendental functions sin u, 
cos u, sinh u, cosh u, to a circle and a rectangular hyperbola with unit semi-diameters, 
with no arbitrary parameter. 


(It is remarkable that these two analogies have escaped the attention of 
mathematicians, possibly from their preoccupation with the cult of general 
mathematical dialectics, during the last three quarters of a century.) 


Simple relations are established 
P, Ps=M, Px = l, 


accents on the suffixes 2’, 4’, denoting that the complementary parameter m’=1—m 
replaces m in the functions thus denoted. 


The substitution p=s-* gives corresponding textbook relations between 
S;, Soy Sy, S,, expressed in terms of snu, cnu, dnu, including Jacobi’s imaginary 
transformation. 


The term “ modulus” is used indiscriminately in textbooks for the real positive 
quantity of the same numerical magnitude as a complex quantity under consideration, 
for the parameters m and k, and for the periods ©, w’. 


In the paper its use is restricted to the first definition, and all quantities with 
implicit complex signs are replaced by their modulus affected by a wholly explicit 
complex sign of the form e, called the generalised sign, 6 being called the phase 
of the complex number. The changes of the generalised sign as the path of 
integration passes the singular points of P~! at the roots 00, 1, m, 0, are established 
in the smallest detail and the results are tabulated against the quantities which they 
affect. In this way the discussion of each case as it arises, often difficult, is avoided. 


Every quantity discussed is shown graphically and extensive use is made of 
the Argand diagram. 


The paper thus far is commended to the attention of students of physics and 
engineering as a logical introduction to the simpler forms of elliptic functions which 
are immediately amenable to numerical applications. The tables of the 
p-functions at the end of the paper avoid the hitherto necessary labour of converting 
to Legendre’s elliptic integrals, or Jacobi’s snu, for which tables are available. 


In Sections 23 to 25 examples are discussed which call for a more extended 
knowledge of the subject. 


In Sections 26 and 27 general methods are used to establish the additional 
theorem and its application to any point in the whole complex plane, with a method 
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of expressing the result in terms of the p-functions and thus making numerical 
applications, by use of the tables at the end of the paper. 


In the last Section it is stated that the zeta and sigma functions may be 
normalised in the same way, and that Lamé’s polynomials in their normalised form 
may be computed and tabulated without prohibitive labour. 


This opens a new field of computation and raises the question of demand 
and supply. In the case of Bessel functions, computation of the simpler forms has 
been followed by more and more generalised forms and now absorbs a predominant 
share of computing activities. It may also be recalled that problems of wind tunnel 
interference have been solved in terms of theta functions, supply in this case 
anticipating demand by more than a century. The question will no doubt receive 
due consideration by computing bodies and institutions. 


| 


